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Abstract. We study the uniqueness and regularity of the steady states of the diffusively 
driven Boltzmann equation in the physically relevant case where the restitution coefficient 
depends on the impact velocity including, in particular, the case of viscoelastic hard-spheres. 
We adopt a strategy which is novel in several aspects, in particular, the study of regularity 
does not requires a priori knowledge of the time-dependent problem. Furthermore, the 
uniqueness result is obtained in the small thermalization regime by studying the so-called 
quasi-elastic limit for the problem. An important new aspect lies in the fact that no entropy 
functional inequality is needed in the limiting process. 



1. Introduction 

1.1. General setting. We investigate in the present paper the properties of the steady 
states of the spatially homogeneous diffusively driven inelastic Boltzmann equation for 
hard spheres interactions and non-constant restitution coefficient. More precisely, we con- 
sider inelastic hard-spheres particles described by their distribution density F = F{v) ^ 0, 
?; G M'^ and we consider the case in which F satisfies the stationary equation 

Qe{F,F) + i^lAF = (1.1) 

for some positive thermalization (or diffusion) coefficient fi > 0. Moreover, assume F has 
a given mass g > and vanishing momentum: 



/ F{v)dv = Q, / F{v)vdv = 0. 



The diffusion operator /iA„F(u) appearing in ( |l.lD represents a constant heat bath which 
models particles uncorrelated random accelerations between collisions. The quadratic col- 
lision operator Qe{F,F) models the interactions of hard-spheres by inelastic binary colli- 
sions where the inelasticity is characterized by the so-called normal restitution coefficient 
e(-) that we shall assume here, in contrast with previous contributions on the subject, to be 
non-constant. This restitution coefficient quantifies the loss of relative normal velocity of a 
pair of colliding particles after the collision with respect to the impact velocity. Namely, if 
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V and denote the velocities of two particles before collision, their respective velocities v' 
and after collision are such that 

(n' • n) = — (n • n) e{\u ■ n\), (1.2) 

where the restitution coefficient e := e(|n • n|) is such that ^ e ^ 1. The unitary vector 
n e §^ determines the impact direction, that is, n stands for the unit vector that points 
from the u-particle center to the t;^ -particle center at the instant of impact. Here above 

U = V — V-),, u' = v' — v'^, 

denote respectively the relative velocity before and after collision. Assuming the granular 
particles to be perfectly smooth hard-spheres of mass m = 1, the velocities after collision 
v' and v'^ are given, in virtue of ( |1.2p and the conservation of momentum, by 

/ 1 + e _ , 1 + e 
V =v — [u-njn, = V-i, -\ — {u-n)n. U-oJ 

The main assumption on e(-) we shall need for our analysis is listed in the following, see 

Assumptions 1.1. Throughout the paper, one assumes the following to hold: 

(1) The mapping r e M+ i-)- e(r) e (0, 1] is absolutely continuous and non-increasing. 

(2) The mapping r e M+ i-)- -deir) ■= r e(r) is strictly increasing. 

(3) There exist a > and 7^0 such that 

e(r) ~ 1 — ar''' as r ~ 0. (1.4) 

The assumption that e(-) is non-increasing can be relaxed and replaced by the more 
general Assumptions 3.1 in [^] (notice that, if e(-) is non-increasing, it is proven in 
Appendix A] that such Assumptions 3.1 are indeed satisfied). In several places in our 
analysis, we shall need slightly stronger assumptions on the restitution coefficient that 
will be properly stated when needed. When no supplementary assumption is specified 



means that the stated result is true under the sole Assumptions |1.1| . Notice that all these 



assumptions will be met by the visco-elastic hard-spheres model which is the most physically 



relevant model for applications [O]. For such a model, the properties of the restitution 



coefficient have been derived in [ |13[ , |24l ]; in particular, e(r) can be defined explicitly by 
the following series 

00 

e(r) = l + ^(-l)Wi r^O (1.5) 

k=l 

where > for any A; e N are parameters depending on the material viscosity. In such 
a case. Assumptions |1.1| are met with 7 = -g and a = ai. In the sequel, it shall be more 



convenient to deal with a second, and equivalent, parametrization of the post-coUisional 



velocities. Fix v and Vi, with u / and let u = u/\u\. Performing in (|L3[) the change 



of unknown a = u — 2{u ■ n)n G provides an alternative parametrization of the unit 



sphere S . In this case, the impact velocity reads \u-n\ = \u\ J and the post-collisional 



velocities v' and v ^ are then given by 

, ^u—\u\a , — \u\a 
v' = v-l3 v'^ = v^ + l3 ^ (1.6) 



1^) . 1] . ™. representation allows us ,o give a precise 



where /3 = = /3 j^|ni ■ 

definition of the Boltzmann collision operator in weak form by 

Q,{fJ){v)^l,{v)dv = 

\ I f{v)f{v,) (^l^iv') + V-K) - - Hv*)) dv* dv (1.7) 

for any test function ip = ip{v). Here, the post-collisional velocities v' and v'^ are defined 
by (|]^. Notice that 

I /i2 I /i2 I i2 I i2 I i2l~^"'7 /-, A I jl — U ■ a 

\v \ + l^^il - \W - \v*\ = -|w| ^ 1 

thus, it follows that (see for details) 

/" Qe{fJ){v)\v\^dv = - [ f{v)f{v.)^e{\u\^)dvdv.^O (1.8) 

where the energy dissipation potential *e is given by 

^3/2 ^1 

*,(r) := — / (1 - e(\/fz)2) dz, Vr > 0. (1.9) 

Notice that, under Assumptions the mapping *e(-) is convex and non-decreasing (see 
again []^]). The functional 

Uf)-= f f{v)f{v.)^e{\u\^)dvdv^ (1.10) 

JM3xM3 

can be seen as an energy dissipation functional for the operator Qe- In particular, multi- 
plying dl.ip by one sees that 

leiF) = QlJLQ 

for any solution F to ( ]l.lp with mass 

Stationary solutions for equation ( pTTj ) in the case of constant restitution coefficient have 
been studied from the mathematical viewpoint by different authors. Existence of such 



solutions was shown in []14|]. The study of moments and tails was described in [11]. 
Uniqueness and stability of these steady states (in the elastic limit) was presented in [20]. 
Different kinds of forcing terms have also been considered in the literature. In particular, 
for the inelastic Boltzmann equation in self-similar variables (corresponding to an anti-drift 
forcing term), stationary solutions correspond to the so-called homogeneous cooling state 
and uniqueness, study of the elastic limit and convergence to self-similarity were presented 
in [p!9[]. Uniqueness of steady states for the Boltzmann equation under the thermalization 
induced by a host medium with a fixed Maxwellian distribution was recently presented in 
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[^] . We mention that the case of dissipative Maxwell molecules has been studied as well 



in nTon and [12]. 



Regarding the existence of stationary states, it is not very difficult to extend the results 



given in [ [140 to non-constant restitution coefficient e(-) and obtain the existence of a steady 
solutions for the diffusively driven Boltzmann equation 

Theorem 1.2. Assume that the restitution coefficient e(-) satisfies Assumption Then, for 
any fi > and any g > 0, there exists a nonnegative F = F{v) e L2(M^) n such that 

Qe{F,F) + iiAF = 



with / F{v) dv = Q and / vF{v) dv = 0. 

Jm.3 Jm.3 



The proof of this theorem follows the path given for constant inelasticity parameter 
in [ |l4l ] and can be deduced from the properties of the solution to the Cauchy problem 
associated to . We refer to Appendix B for the main steps of the proof of the Theorem 



1.2. 



1.2. Scaling argument and formal limit A — ^ 0. Let us discuss the main concern of the 
present work, namely, proving the uniqueness of solutions to ( |l.ip 

Qe{F,F) + fiAF = 

in the weak thermalization regime, i.e. when the diffusion parameter fi is sufficiently small. 
In order to understand this regime and the strategy fix ^ > and denote by F a solution 
to ( |1 . ip with given mass g and vanishing momentum. Introduce the following rescaled 
solution 

Gxiv) = X^F{Xv), A>0 (1.11) 
and the rescaled restitution coefficient 

= e(Ar) Vr > 0. 

Since 

X^Qe{F,F){Xv) = Qe,{Gx,Gx){v) and X^A,F{Xv) = A,Gx{v) 
for any v G R^, one gets that Gx is a solution to the rescaled stationary problem 

Qe,{Gx,Gx) = -^A,Gx. (1.12) 

In other words, for any A > the rescaled distribution Gx is a solution to the steady dif- 
fusively driven Boltzmann equation with thermalization coefficient /i/A^ and restitution 
coefficient ex (notice that ex still satisfies Assumptions |1.1[ ) . For any A > 0, the solution 
to dl.lD is unique if and only if the solution to ( ]1.12| ) is unique. Such a scaling is partic- 
ularly interesting because, in addition to preserve mass and momentum, when A — ^ the 
rescaled restitution coefficient ex{r) converges pointwise to the elastic restitution coeffi- 
cient lim^^o ^x{f) = 1 for any r ^ 0. Consequently, one formally expects that 

Qe,(/,/)^Qi(/,/) as A^O 
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where Qi (/, /) denotes the classical Boltzmann operator for elastic interactions. This 
means that the dissipation of energy is expected to vanish as A — ^ 0. Formally, one sees 
that if > is kept fixed the right side of ( ]1.12| ) will be infinite in the limit A — > 0. In other 
words, the thermalization fi has to be reduced to compensate the loss of dissipation, i.e. 
one has to choose a diffusion coefficient fi = /ix depending on A such that liniA^o A'a = 0. 
Intuitively, it make sense to look for a parameter that keeps the solution's energy 

£x = - f Gx{v)\v\Uv 

of order one in the limit A — > 0. Let us investigate the correct scaling by multiplying ( ]1.12| ) 
by |wp and integrating over M.^ to get 

6nxg = 1^e^{Gx) 

where I^x is the energy dissipation functional associated to the rescaled restitution coeffi- 
cient ex given by 

Ie,{f)= f f{v)f{v,)^,{\M'')'^vdv,. (1.13) 

Since is convex, a simple use of Jensen's inequality yields QfixQ ^ ^?^*e (-^^"^a) • More- 
over, 

3+7 

cir 2 

^e{r) ^ as r ~ 0, (1.14) 

4 + 7 

2 

from which it follows that \^8x = 0(/U^^^) as A ~ 0. Consequently, to keep the kinetic 
energy £x of unit order we must have 

^i■.= ^Jix = A'+^. 

For such a scaling, equation ( |1.12p becomes 

Qe,{Gx,Gx) + \^^vGx = Q. (1.15) 
Note that with our choice of ^x the limit Gq as A — ^ of Ga, if it exists, has to satisfy 

Qi(Go,Go) = 0. 

In other words. Go is a suitable Maxwellian with same mass and momentum that Ga. 
Moreover, using the dissipation functional 

6f? = 4/ Gx{v)Gx{v.)^e{\^\v-v^\'')dvdv, (1.16) 

one expects that the limit Go satisfies 

(iQ= I Go{v)Go{vM\v-vA'')dvdv^ (1.17) 

with Co(?'^) = limA^o ^*e(A^r^) (several properties of such energy dissipation functional 
are investaged in Appendix A). With this observation, it is not difficult to prove that the 
unique possible limit as A — ^ of Ga is the Maxwellian distribution 
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for some explicit temperature Q determined by the above identity (|1.17[) 



The limiting Maxwellian M{v) is called the quasi-elastic limit for this problem. With 
this knowledge we will prove uniqueness of solutions for the problem ( |1.15| ) when the 
rescaling parameter A is small (lying in an explicit interval) . In this sense our uniqueness 
result will be valid in the weak thermalization regime. 

1.3. Main results and strategy. Let us state precisely the main problems we wish to ad- 
dress in this document: 

(1) Prove that any solution G\ to ( |1.15P satisfies 

lim Gx=M 

in some suitable sense. Specifically, find a suitable Banach space X such that 
G\£ X for any A > and limA^o ||Ga - M\\x = 0. 

(2) Prove that solution G\ to ( |1.15[ ) is unique, at least in the weak elastic regime. That 



is, determine A^ G (0, 1) such that &\ reduces to a singleton as soon as A G (0, A^]. 
We use the symbol ©a to denote the set of solution G\ to ( p..l5p with given mass 
and vanishing momentum. 
(3) Provide quantitative answers to the two previous questions. More precisely, find 
the rate of convergence of G\ towards M as well as some estimate for the param- 
eter A^. 



The first question is answered with the following theorem (see Theorem |4.1| for a de- 
tailed statement) that can be interpreted as a quasi-elastic limit result. 

Theorem 1.3. If e{-) is of class C" with m ^ 3 (with some additional regularity properties), 
one has 

lim ||Ga - A^IIh* = VA; ^ 0, V^G[0,m-2]. 

A — >0 ^ 



The limit M. is the Maxwellian given by ( |1.18| ) with an explicit temperature Q given by ( |4.1p . 



The convergence also holds in exponential weighted L^-spaces. 

The proof of the above result is based upon a compactness argument and requires a 



careful investigation of the regularity properties of the solution to (1.15). Our approach 



for the study of regularity of solutions to (TTTsJ) differs from the related contributions on 



the matter [ ]20| , [T^ ] where the regularity of steady solutions is deduced from the properties 
of the time-dependent problem (namely on the propagation of regularity combined with 
the damping in time of the singularities for solution to the time-dependent problem, see 
^]). In contrast with these results, our methodology is direct and relies only on the steady 



equation ( 1.15p . It is not difficult to prove by a bootstrap argument that any solution 
Gx to ( |1.15 ) is smooth, however, it is more delicate to obtain regularity estimates which 
are uniform with respect to the parameter A > since the diffusive heating in ( |1.15P is 
vanishing in the limit A — ;> 0. On the basis of new regularity estimates of the collision 



operator (see Theorem 2.5), we can prove the following proposition. 



Proposition 1.4. Under the regularity assumptions on e(-) of Theorem |j.3| one has 

sup ||Ga||h* < oo V/c ^ 0, £ £ {0,m — 1]. 
Ae(o,i] 



Theorem |1.3| serves as a fundamental brick to prove the main result of the paper. 

Theorem 1.5. Under suitable regularity assumptions on e(-) there exists A^^ G (0, 1] such that 
the set &x of solutions to ( |1.15[ ) with given mass g and vanishing momentum reduces to a 
singleton for any A G [0, A''']. 



Theorem |1.5| can be interpreted as an uniqueness result in the quasi-elastic regime where 
A is small. This theorem, however, can also be interpreted as a weafc thermalization unique- 
ness result since in this regime the diffusion parameter /i is small as well. 

Theorem 1.6. Under suitable regularity assumptions on e(-), there exists /i^^ > such that, 
for any fi G (0, fi^ the steady problem 

Qe{F,F) + iiAF = 

admits an unique solution F for a given mass and vanishing momentum. 



The proof of Theorem |1.5| follows the strategy of [ ]20| | (see also [ ]T^ ] and [^). Essen 



tially, it is based on the knowledge of the quasi-elastic limit problem and on quantitative 
estimates of the difference between solutions to the original problem and the equilibrium 
state as A — )> 0. More precisely, let us consider two steady solutions G\,Fx G &\. Set then 
Hx = Fx — Gx and define the linearized elastic Boltzmann operator around the limiting 
Maxwellian M 

^i{h) = QiiM,h) + Qi{h,M). (1.19) 
Observing that is a S3mimetric operator one recognizes 

^i{Hx) = (qi{Hx,M) - Qe,{Hx,M)] + (qi{M,Hx) - Qe,{M,Hx) 



+ Qe.iM - Fx,Hx) + QeAHx,M -Gx)]- X<AH, 



where we used that Qe^{Fx, Fx) — Qex{Gx,Gx) = X'^AHx- Assume that there exist two 
Banach spaces X and y independent of A such that 

\\QeAL9)\\;, + \\Qe,{9j)\\;,^Gi\\f\\y\\g\\y, (1.20) 

and 

||Qi(/,A^) - Qe,(/,A^)||;t- + WQiiMJ) - QeAM,f)\\;^ ^ C2X^f\\y, (1.21) 
and also, 

WAHxWx < G-iWHxWy VA > (1.22) 
for some constants Ci, C2, C3 > 0, p G (0,7) independent of A. Notice that ( |1.22p is too 
restrictive for a general function / G 3^, it is only assumed for any difference Hx- Then, 

\\^i{Hx)\\x ^ Ci\\Hx\\y{\\M - Gx\\y + \\M - Fx\\y) + {G2 + G3)X^\\Hx\\y. 
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If M is the universal limit of the family &x, if one is able to prove that 

lim - Gx\\y + \\M - FxWy) = (1-23) 

then, for any e > there exists Aq > such that 

\\^i{Hx)\\x ^e\\Hx\\y. 

The strategy ends by proving that there exists a subspace y C y containing the net 

{^^a}ag(o,i] such that 

\\^iih)\\x co\\h\\y yh€y. (1.24) 

Thus, for any e > there exists Aq > such that 

coll-ffAb ^ell^Ab VAg(0,Ao). 
This proves that Hx = for any A e (0, Aq). Notice that Aq would become explicit if we are 



able to make explicit the rate of convergence of (|1.23|) . 



To summarize, the proof reduces to find Banach spaces X and y for which the above 
equations ( ir20| )-(|L24l) hold. We warn the reader here that ( |1.24[ ) will have to be slightly 



modified because a priori the energy of the difference Hx is not necessarily zero. This 
technical detail is overcame by introducing a suitable lifting operator of ^i, see [ ^0[ ] for 
the original implementation of this idea. We can already anticipate that the strategy will 
be applied to the following weighted L^-spaces 



X = L^inia) = L^{R^; ma{v)dv) and y = L\{ma) = L^{R^; ma{v) 1 + \v\^ dv) 
where the exponential weight function nia is given by 

ma{v) = exp {a\v\) a ^ 0. 

The most technical parts of the proof will be the A-uniform regularity of Gx and the 
continuity estimate ( Il.21| ) with respect to the restitution coefficient. These aspects have 



been proved for constant restitution coejficient in []20P, however, their extension to the case 
of a variable restitution coefficient will be delicate and require a series of new technical 
results. 



Finally, explicit estimates on the rate of convergence of the rescaled solution Gx towards 
the elastic limit Ai can be found a posteriori by seeking a nonlinear inequality satisfied 
by II Ga — A^lly- More precisely, we shall prove that there exist some explicit constants 
Ci,C2 > such that 

IIGa - A^b ^ CiXP + C2\\Gx -Mfy VA e (0, 1]. 

Combining this estimate with the convergence of Gx towards Ai will lead to the existence 
of some A^" G (0, 1] such that 

WGx-MWy^GsXP VAe(0,At], 
for a suitable exponent p > and explicit constant C3. 
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We stress here that in contrast with the reference []20|], the present manuscript provides 
an approach which does not rely on entropy estimates. Consequently, it does not require 
neither exponential pointwise lower bounds nor strong regularity properties in the steady 
state. In particular, it is well-suited for problems in which no regularity of the steady 
solution is available, see for an example of this situation. 



1.4. Organization of the paper. The plan of the paper is the following. In Section 2, we 
establish new regularity estimates of the collision operator Qf generalizing known results 
for the elastic case [ |T6| , |27| , |22[ ] and for the inelastic case with constant restitution coeffi- 
cient [ ]T^ , [T7| ] . Section 3 is devoted to study regularity properties of the steady solution 
Gx G ©A- In particular, we study moments [|ll], ^ ^, general weighted Sobolev regu- 
larity given in Proposition |1.4| and a technical result on the difference of two solutions, 
see Proposition 3^. Moreover, we also address in this section the fundamental problem 
of the continuity properties of Qe^ with respect to the inelasticity parameter, proving the 
convergence of Qexif^d) towards Qi{f,g) as A in different norms. Several of these 
results are non-trivial extensions of those in [ ]T9| ] given for constant restitution coefficient. 
Others, like Proposition |3.1Cl| are new. Section 4 contains the main results of the paper. 



namely the elastic limit result Theorem |1.3| , the uniqueness result Theorem |1.5| and their 
quantitative versions. In Appendix A, we present several technical results used throughout 
the paper and Appendix B contains a proof of Theorem |1.2| on existence of steady profiles 
adapted from [Q. 



1.5. Notation. Let us introduce the notations we shall use in the sequel. Throughout the 
paper we shall use the notation (•) = y^l + | • p. We denote, for any p G [1, +oo), G M 
and weight function txx : — ^ M+, the weighted Lebesgue space 

= |/ : m3 ^ R measurable; := \f{vW{vrzu{v)dvy'' < +oo| . 

Similarly, we define the weighted Sobolev space W^5'^(tz7), with /c e N, using the norm 

\\s\sik 

The symbol denotes the partial derivative associated with the multi- index s G N^: 9** = 
dyj^d^'^d^^. The order of the multi- index being defined as |s| = si + 82 + 53- In the particular 

case p = 2 we denote H,^'(ti7) = wJ5'^(ti7) and whenever w{v) = 1, we shall simply use H.^. 
This definition can be extended to for any s > by using Fourier transform. 




2. Regularity properties of the collision operator 



The smoothing properties of the gain operator Q^{f,f) have been investigated in our 
previous contribution \^]. However, for the results we have in mind, we shall need the 
regularity of the bilinear operator Qfif, g) rather than the one of the quadratic one. To do 
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SO, we shall use a slightly different approach than the one we used in . In particular, our 
main purpose here is to extend []6[ Theorem 4.1] to smooth kernel that are not compactly 
supported: in such a case, the price to pay for the control of large velocities consists in 
additional moments estimates. Precisely, using the notations of Appendix A, let B{u, a) be 
a collision kernel of the form 

B{u,a) = <i>{\u\)b{u ■ a) 

where <!>(•) ^ and b{-) ^ satisfies (A.3) and (A.4) of the Appendix. Then, one can 
define the following operator Fb by 

Tb{'p){x) = / B{z+ae{r)uj,ae{r))(p{ae{r)u)+z) d7r{z), x = rw, r ^ 0, G S^, (2.1) 

where dvr is the Lebesgue measure over the hyperplane perpendicular to oj and ae(-) is 
the inverse of the mapping s s/3e(s) while the kernel B{-, •) is given by 



with '&e{-) defined in Assumption |1.1| (2) and denoting its derivative. The operator 
can be seen as an inelastic version of the so-called cold thermostat operator investigated 
in \^ (and originally derived in the seminal paper [Jilt]) and plays a crucial role in the 
smoothing properties of the gain operator Q+ because of the representation formula 



QheU,9){v)= f{z)[{t,oTBot,)g]{v)dz \/f,g (2.3) 
where [t^'ip]{x) = ^{v — x) for any i;, x G and test-function (see for the derivation 



of ( |2.3p ). Before investigating the regularity of the full gain operator, we shall first deal 



with that of the cold thermostat. 

2.1. Regularity properties for cut-off collision kernels. For this section we assume that 
the kernel B{u,a) satisfies: 

GC-(0,oo), 6(.)GCo°°(-l,l) and m = {ltlltle, ^^'^^ 
for some e > 0. We introduce the following definition: 



Definition 2.1. We shall say that a restitution coefficient e(-) satisfying Assumptions |TT 
belonging to the class Km for some integer m ^ 1 if e{-) G C™(0, oo) and 



IS 



sup re^''\r) < oo yk = l,...,m (2.5) 
where e^''\-) denotes the k-th order derivative of e{-). 

Remark 2.2. For the physically relevant case of visco-elastic hard-spheres, the restitution 



coefficient e(-) is given by ( |1.5p but admits also the folloM^ing implicit representation (see 

i 3 

e(r) + ar5 e5 (r) = 1 Vr > 



11 

for some a > 0. Then, it is possible to deduce from such representation that e(-) belongs to the 
class Km for any integer m ^ 1. 

Under these assumptions we have the following generalization of Lemma 4.6]. 

Lemma 2.3. Assume that e(-) belongs to the class Em for some integer m ^ 2 and that the 
collision kernel B{u, a) satisfies assumption ( |2.4| ). Then, for any ^ s ^ m — 2, there exists 
C = C{s, e, e) such that 

l|rB(/)|lH^+i^C||/llH^ Vr?^0 (2.6) 

with fi{s) = s + 4 and where the constant C{s, e, e) depends only on s, on the collision kernel 
B and the restitution coefficient e(-). 

Proof There is no loss of generality in assuming that s is an integer. The proof is divided 
into five steps. 

• First step: change of variables Recalling Lemma 4.6], we first define 

f^(/)(ru;) = TBif)ia-\r),co) = rB(/)(r/3e(r), a;) (2.7) 

so that 

TB{f){ruj) = / B{z + ruj,r)ip{ru + z) dTT{z). 



We begin proving the result for Tb instead of T^, that is 



^Bif) ^ C{s, B, e) II/IIh. ^ , Vr? ^ (2.8) 



The proof of this estimate follows the approach given in p2| , Theorem 3.1] where a similar 
estimate has been obtained, for /i(s) = 0, under the additional assumption that <!>(•) has 
support in [e,M] with M < oo. Our proof will consists essentially in proving that the 
weighted estimate (i.e. with fi = fi{s) > 0) allows to take into account large velocities. 

• Second step: Estimates on the radial derivative o/r^. We introduce the radial Fourier 
transform TZT and the Fourier transform in with the formulas 

TZT [f] (gw) = (27r)-i/2 / exp(i^>r)/(r«;) dr , T [f] {0 = (27r)-3/2 f exp{ivOf{v) dv 
and, for any measurable mapping g, we define the H^+^(S^ x M) norm of g as 



|2 



II^IIh^+1(s2xr) 

Then we compute, 
TIT 



S2 



{r)^TBf (qw) = {2tt) 2 I exp{iQu ■ w){u ■ w)^B{u,\u ■ w\)f{u) du 

= 27r^[/(-)G^(-)] (gw) 
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where Gw{u) = B{u, \u ■ w\){u ■ w)'^ for any u G "M?. Then, setting ^ = gw, since Agdw 
1^1^^ we get 



H^+l(§2xI 



2tt 



/(•)Gi(- 



(6 



Now, with fi = fx{s) = s + 4, introducing g{v) = fiv){v)'^ and Gw{z) = {z)^^^Gw{z) we 
can write the above as 

2 



H^+i(S2xI 







5(-)G^e (•) 


(0 


/R3 




1«T . 





de. 



Using [ ]2^ , Lemma A.5] we can replace Q24 Eq. (3.7)] by 



sup 



B{z, \z • wl 



{z ■ wY 



(2.9) 



• Third step: Control of the angular derivatives. In order to adapt the analysis of []22[], it 
suffices to check that there are two positive constants a,b > such that 

Supp(;B) C [a, 00) X [6,00). 

We already saw this is the case with our assumption on $(•) and We can straightfor- 
wardly apply the reasoning of the op. cit. to get that the j-th angular derivative of TbJ 
can be estimated by the radial derivative of r^^ / where the new kernel Bj is given by 
Bj{z, g) = B{z, g)zj. This finally leads to (2.8) with 

.{z-w)'^ 



C{s,B,e) 



Cs{a,b) sup sup 

|j^|sSs+l 



B{z, \z ■ w\)z'' 



(2.10) 



Fourth step: Let us check that the above quantity is indeed finite, i.e. 

.{z-w)'^^ 



Cs{B,u) 



sup 



B{z, \z ■ w\)z'' 



{z)v^ 



< 00 



Hs+2 



for any multi- index u with ^ s + 1. Observe that, because of our cut-off assumptions 
( |2.4| ) together with the fact that 6(1 — x) = for small values of x, the kernel B{z, g) 
vanishes for small values of \z\ and g. Thus, for a given u with ^ s + 1, it suffices to 
investigate the regularity properties of the above mapping 

(z ■ w)^ 
: z^B{z,\z-w\)z^ ^, ^ ' 

for large value of z (uniformly with respect to w). From the definition of Cs(^,z/), one 
needs to compute s + 2 derivatives of which explains the restriction s ^ ?ti — 2. Recall 
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that, for \z\ > e, the expressions of B{-, •) and $(•) yield the following 

8b{l - 2\z ■ w\^) z" {z ■ w)'^ 
"^^^ " niz-wl) 

where z = z/\z\, 

nir) =r(3^{r){l + i9'^{r)), r ^ 0. 

We recall that 6(1 — 2x^) = for |x| ^ 6 for some 6 > 0. In particular, for \z\ > e, 
Fw{z) / =^ r := \z ■ w\ > 6e. Since /3e(-) G (1/2, 1] and ^ 0, it is easy to check that 

Since /i — s = 4, this proves that sup^g§2 ||-^w(-z)||l2(k3) < oo. One proceeds in the same way 
with the z-derivatives of Fy^{z). It is clear that any z-derivative of the rational expression 

z''{z-w)'^ 



Rw{z) :-- 



has a faster decay (for |z| — ;> oo) than Rw{z)- Therefore, the crucial point is the control of 
the derivatives of nj^. It turns out that 

n{r} 

n'ir) _l 2/3^(r) ^'^{r) _1 2e'(r) 2e'{r) + re"{r) 
nir) /3e(r) ^ l + t?'g(r) ~ r ^ 1 + e(r) ^ l + t?'g(r) ' 



Now, our assumption ( |2.5p on the restitution coefficient e(-) implies easily that V.' /V. G 

L°°([(5e, oo)) and, as a direct consequence, 

^^eL~([&.oo)). 

Similar calculations show that, for any k = l,...,m, -^yjy^ G L°°{[6e, oo)). Tedious 
but simple calculations show then that any z-derivative of F^(z) can be controlled by 
l/l^j/i-kl+i for large \z\. This is enough to prove that Cs{B, v) < oo. 

• Final step: turning back to the original variables. Following [^, it remains now to 
deduce estimates on Tsf from Tsf which are linked by formula ( |2.7[ ). Using polar coor- 
dinates 

/■oo 

l|rB(/)||L = E / I'Me'ief'do / \diTB{f){Q,u^)\'duj 

where, see for details, one can check that for any \j\ ^ k the function Fj{g) can be 
written as 

FjiQ) = pMHq), . . . ,#(^'))(1 + ^«(^>))-"^. (2.11) 
Here Pj {yi, ... ,yj) is a suitable polynomial, nj G N and i^i^-* denotes the p-th derivative of 
i?e(-)- Because of our assumption on e(-) (more precisely, because limsup^_^3o ^e\Q) < oo), 
we see that sup^gj-g^oo) = Cj < oo for any \j\ ^ k. Thus 

l|rB(/)||H^ ^c,||f^(/)||H^ (2.12) 



14 R. J. ALONSO & B. LODS 

where is an explicit constant involving the L°° norm of the first s-th order derivatives 
ofae(-)- □ 



Proposition 2.4. Let B{u,a) = ^{\u\)b{u ■ a) be a collision kernel satisfying ( |2.4{ ) and e(-) 
be in the class {m ^ 2). Then, for any ^ s ^ m — 2, 

with constant C{s, B, e) and /i(s) = s + 4. 



Proo/ One uses the representation formula ( ]2.3p together with Minkowski's inequality to 
get that 

l|2i3e(/'fi')llH''+^ ^ / 1/(^)1 11*2 °rot^5r||jj«+idz. 

Now, since H^zV'IIh^ ^ (•2)'^||V'IIh^ for ^^Y ^ ^ for all G N and any k ^ 0, one easily 

fc fc 



deduces from Lemma 2.3 the conclusion. □ 



2.2. Regularity properties for hard-spheres collision kernel. We now use the previous 
result for smooth collision kernels to estimate the regularity properties of Q^{f, g) for true 
hard-spheres interactions. We shall combine Theorem A. |l| of the Appendix together with 
the estimates of the previous section to get the following: 

Theorem 2.5. Assume that e(-) belongs to the class with m ^ 2. Then, for any e > and 
any r] ^ 0, there exists Ce = C{e, e, t]) such that 

\\Qt{f,9)\\rf-+^ ^Ce\\g\\H^ ll/Li IbllH- 



+ e IblLi , II^VIIl^ + ll/ILi \\d'g\\L2 m =s + l^m-l. (2.13) 

\ ij+i 17+1 f+i f+i/ 

Proof Notice that, for hard-spheres interactions, one has B{u,a) = ^{\u\)b{u ■ a) with 
<l>(|n|) = \u\ G L°^-^^ and b{s) = I/Att for any s € (—1, 1). In particular, for any rj ^ 0, both 
the constant C2,i,,,,i(6) and C2,2,r),i(^) appearing in ( |A.5P are finite. Let us now Bx -q ^ 
and e > and split the kernel into four pieces 

B{\u\,u-a) = <^s{\u\)bs{u ■ a) + <^si\u\)bR{u ■ a) 

+ ^R{\u\)bs(,u ■ a) + <^R{\u\)bR{u ■ a) (2.14) 

with the following properties: 

(i) bs and are smooth satisfying the assumptions of the previous section. 

(ii) bR{s) := ^ — bs{s) is the angular remainder satisfying 

C2,i,-n,i{bR) ^ e and C2,2,,,,i(6ij) ^ e. 

(Hi) ^r{\u\) = \u\ — ^>5'(|?i|) is the magnitude remainder satisfying 

e 



iC2,l,rj,l{bs)+C2,2,v,l{bs))' 
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Notice that, in contrast to previous approaches, the last point is made possible because 
= 1^1 for large \u\ which makes <i>/j compactly supported. Thus, on the basis of 
relation ( |2.14D , one splits Qf into the following four parts using obvious notations. 

We shall then deal separately with each of these parts. First, we know that 

\\Qssif^9)\\uf,+' ^ C'„,„,e||9l|H^,+^||/|lLi^_^^ Vm,n. 
Second, let us estimate Q^^. Since 

for any multi-index i with \i\ ^ s + 1, one gets 

We treat differently the cases \£\ = s + I and \£\ < s + 1. According to Theorem A. |l| if 

1^1 ^ s one has for any \u\ ^ \£\ 

rj+l 77 + 1 

where we used the assumption (ii) with the fact that ||$5||j;^oo^ ^ 1. Recall the general 
estimate 

WoWl^ ^ ^e\\g\\L^ yk^o, ye>o (2.15) 

where the universal constant tq is given by = ||(-) 2 ||^2 < 00. Taking for simplicity 
9 = 1/2 and since |£| ^ s, 

|£|<s+li'=0 ^ ^ 

for some constant > depending only on s. In the case |^| = s + 1, argue in the same 
way to obtain 

for any < < |£|. If = one still has 

additionally, for = £ we use Theorem A. [l| with {p, q) = (2, 1) to get 

\\QUd'f,9)\\L^, ^ C2,2,vAibR)\\9\\Ll^, Wd'fh^^,- 



s + 1. 
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Therefore, 

Third, argue in the same way using the smallness assumption (ii) to prove that 

\\QRRif^9)\\^^+- ^ 

Finally, the estimate for follows from the fact that ||<I>ij||L°° is small, 
II2^^5(/'5)IIh^+^ ^ (II/IIh^,, IbllHj; + llffllLi PVIIl^ + II/IIli Wd'gh^^ 



\f\e\ = s + i. 



Combining all these estimates and replacing AsE to e we get (2.13). 



□ 



Remark 2.6. Recall that, by virtue of our scaling argument, we will have to apply the above 
regularity result for the scaled restitution coefficient ex. Arguing as in [^, Corollary 4.14] 
we can prove without major difficulty that sup;^g(o,i] < oo where Ce^ is the constant 
appearing in (|2.13p for the scaled restitution coefficient e\. 



3. Properties of the steady state 
The purpose of this Section is to establish all the general a posteriori properties of the 



family {G\)\ of solutions to ( |1.15p that will be necessary to establish the uniqueness result. 
Of course, this analysis will require fine properties of the collision operator Qe^ associated 
to the rescaled restitution coefficien t ex, in particular, its behavior as A — ^ 0. In all this 
section, Gx denotes any solution to ( |1 . 1 5p with A e [0, 1]. There is no loss in generality in 
assuming from now on that 

g= [ Gxiv)dv = 1 VA G (0,1]. 

We shall define, for any A G (0, 1] the solution set: 

<3x = ^Gx £ lI;Gx solution to ( |1.15P with 

/ GA(w)dw = land / vGxiv)dv = o\. (3.1) 
Recall that our choice of scaling implies that for any Gx G &x, the energy identity is given 



by, see ( |1.16[ ) 



^ = T^I Gx{v)Gx{v,)^e{Av-vA'')dvdv, VAg(0,1] (3.2) 
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where *e has been defined in (L9). We deduce from ( 1.14 ) that, for fixed r > 0, 

^ ^^(A^r^) ~ — ^r^+T as A ~ 0. 



A3+7 ' 4 + 7 

Intuitively, one gets then that, for A ~ 0, 

6^-^ [ Gx{v)Gx{v.)\v-v.\''+''dvdv^. 

4 + 7 

The use of Jensen's inequahty proves that the moment of order 3+7 of Gx remains bounded 
uniformly with respect to A 

ms+y{X) = 0(1) 

2 

where the moments are defined as 

mp(A) = / Gx{v)\v\^Pdv p^l. (3.3) 

Existence of higher moments for Gx is the objective of the following section, see also 
Lemma A. ^ in the Appendix, which properly justify above computations. 

3.1. Moment estimates. Recall that our choice of scaling is such that 

sup mi (A) = sup £x = Smax < 00. 

OsgAsSl O^A^l 

By a simple induction argument, this actually implies that all the moments of Gx are 
uniformly bounded. 

Proposition 3.1. For any p ^ 0, there exists Cp > such that 

sup mp(A) ^ Cp. 

O^A^l 

Proof. Let p ^ 1 be fixed. Multiplying Eq. ( |1.15p by ^'(w) = \v\'^'p and integrating over 
we get 

-A^ / Gx{v)A\v\^Pdv= [ Qe,iGx,Gx)iv)\v\^Pdv. 



Since A|i'p*' = 2p{2p + l)|upP ^, using Lemma B. [T], there are two positive constants 
kp,Ap>0 independent of A such that 

-2p(2p + l)A'^mp_i(A) ^ -/cp^.m 1 (A) + Apmi (A) mp(A) VA ^ 0. 

2 

Since mi (A) ^ \/£x ^ V^ma^ for any A G (0,1), we see that there are two positive 

2 

constants Ci,p, C2,p > independent of A such that 

m ^1 (A) ^ Ci,pmp(A) + C2,pmp_i(A) VA G (0, 1], Vp 1. 

Both sup;^g(o,i) ™i(-^) = "^max ^ud sup;^g(-Q i] mo(A) = 1 are finite, thus, a simple induction 
yields the conclusion for any p G N. The result extends then to any parameter p ^ by 
interpolation. □ 
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Proposition 3.2. There exist £^min > and cq > such that inf;^£[o ij £\ = <?min > and 
Gx{v^)\v-v^\dv^^co{v) V^;GM^ VAg[0,1]. 

Proof. Note that ^e{x) — Cx^/^ as x — ;> oo for some positive C > 0. Additionally, using 
( |1.14[ ), there exists a positive constant K > such that *e(?'^) ^ Kr^^'^ for any r > 0. 
According to (]3.2p, it follows 



6€K 



Gx{v)Gx{v,)\v-v,f+^dvdv,. 



Therefore, 



inf m3+^(A) 
A>o V 



c> 0. 



(3.4) 



Knowing ( |3.4[ ), it is a standard procedure to deduce the result from Proposition [3.1| . □ 

Proposition 3.3. There exist positive constants ^ > and M > such that any solution G\ 
to dl.lSp , with A G (0, 1], satisfies 

[ Gx{v)exp(A\v\'^)dvi:M. (3.5) 



Proof The proof follows the lines of the analogous result [ ]Tl| , Theorem 1] for constant 
restitution coefficient. It consists in proving that there exist K > such that 



sup mp(A) ^r( ^p+i ) KP Vp^ 1 

AG(0,1] 



(3.6) 



where r(-) is the gamma function while mp(A) is defined in ( |3.3| ). In order to prove ( |3.6p 
note that, 

-2p(2p+l)A^mp_i(A) = / Qe,(GA,GA)(t;)bP^d7; Vp ^ 1, A G (0, 1]. 
One can estimate the right side thanks to Proposition 2.7], 

/ Qe,(GA,GA)(t;)b|2Pdt;^ -(l-Kp)m i(A) + /tp5p(A) 



where 



'5p(A) = ^ ( ^ ) (mfe+i/2(A) mp_fc(A) + mfc(A) mp_fc+i/2(A)) 
Here [^^] denotes the integer part of and k^G (0, 1) is independent of A and satisfies 



0(l/p) as p 



oo. 



P. 

. Then, one sees that [ ]Tl| , Equations (4.6) and (4.11)] hold with 



= A"^ G (0, 1]. At this point, we can resume exactly the proof of [ ]ll| ] noticing that all 
the estimates there are uniform with respect to the coefficient fi appearing in front of the 



thermal bath. In other words, we obtain ( 3^ with a positive constant K > which is 
independent of A. This is enough to get ( |3.5p . □ 
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One actually can make more precise the above estimates by evaluating the difference 
of two solutions to ( |1.15p . A simple adaptation of [ l20[ Proposition 2.7, Step 1] gives the 
following estimate. 

Proposition 3.4. For any s G [0, |] there exist some positive constants > and Ms > 
such that 

[ \Gxiv)-Fx{v)\expirs\vndv^M4Gx-Fx\\L] VA G (0, 1] (3.7) 
for any Fx, Gx G ©a- 

3.2. Sobolev estimates. We prove now that the family {Gx)x is uniformly bounded in any 
Sobolev norms H^. We begin showing uniform L|-estimates of Gx for sufficiently small A. 

Proposition 3.5. For any k ^ 0, one has := sup;^g(Q ^ ||Ga||2,2 < 



OO. 



,2k 



Proof. First, observe that for any test function \Ij{v) integration by parts yields 

AGxiv)Gxiv)ijiv)dv = [ \VGxiv)\'^i^iv) dv - I- [ Gx{vfA'iP{v)dv. 

Fix k ^ and multiply equation ( |1.15p by Gx{v){v)^''. Apply above identity to ip{v) = { 
and use the inequality A'ip{v) ^ 2k{2k + l){v)'^^'^^^'^ to obtain 

A^llVGAlli^ ^ / Qe,{Gx,Gx){v)Gx{v){vf''dv + i2k^ + k)X^\\Gx\\l2 . (3.8) 

Applying Corollary 4.14] with p = 2 and t] = k,we see that there exist G (0, 1), z > 
and Ce^ > depending only on ex such that for any 6 > 0, 

I Qt{Gx,Gx){v)Gx{v){v)^'' dv ^ Ce^^^-^IGAll^r^llGAllJ^i'"') +<5||GA||ii IIGaII^. . 

Jigs 2+fe ^k+1/2 

Same reasoning as in Corollary 4.15] shows thatQ 

sup Ce^ < OO. 

Ae{o,i) 

Therefore, there exist 9 G (0, 1) such that for any A G (0, 1) and 6 > one can find some 
Ks > independent of A for which it holds 



Qt{Gx,Gx){v)Gx{v){vf' dv 



^Ks\\Gx\\lt''\\Gx\fJ^y^ + S\\Gx\\^JGx\\l. ^ VA G (0, 1) 



"^With the notation of Corollary 4.15], one can prove that for any compact interval / C (0, oo) it 
ows that maxfc^o.i 

liniA^o Ce, = Co > 0. 



followrs that maxfc^o.i \\D'' Ge^{-)\\L^ (i) = 0(1) as A ~ where Ge^(r-) = — — In particular, 

(1 + t'exC^'jjPe;, yr) 
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Second, estimate the loss term thanks to Proposition ^j^. Indeed, 

Q-^{Gx,Gx){v)Gx{v){v)^'' dv = I \v - vAGI{v)Gx{v.){v)^'' dvdv, 

^co f Gl{v){vf>^+Uv = c,\\G^\\l, 



k+l/2 



Thus, plugging the previous two estimates into ( |3.8[ ) 

A^llVGAlli^ ^ KsWGy^W^^.^'WGxf^:^^ + ('^IICaIIli -co) IIGa"^ 



2+fe "/ " ""■^fe+1/2 



+ {2^ + k)X'\\Gx\\l2 VAg(0,1]. 



Using the notation of Proposition p . 1| and choosing 5 = co/2C2+fe, one sees that there exists 

Ck = Ks&f^'^'^ > such that 

>^''\NGx\\l2 ^CkWGxf^t'^ -^WGxWl^ +i2e + k)X"'\\Gx\\l2 VAg(0,1]. 

k Z fc + 1/2 

In particular, 

'-^WGxWl. ^Ck\\Gxft'^ + {2k^ + k)X'\\Gx\\l. . 

2 ■^fe + l/2 ^k-l 

The case k = follows directly from this estimate, i.e. 

sup ||Ga||l2 < oo. (3.9) 

AG{0,1] 

Assume now A; ^ 1. For any > 0, it can be checked that 

||Ga||^2 ^ R^'^'^WGxWl^ + R-HGxWI^ . 



k+l/2 



1 /3 

Hence, choosing R = (^^X"'{2k'^ + k)^ we get 

-|^||GA||r2 ^ CfellGAll^'2" + Bk{X)\\Gx\\'j^2 
4 -^k+1/2 -^k ^ 

with Bk{X) = {2k'^ + k)X'R^^~'^. In particular, using ( |3^ ) there exists some positive 
constant := -Bfc(l) supAg(o,i) ||Ga|||2 > 0, independent of A, such that 

^||GA||i2 ^Ak + Cu\\Gxt^2''^ VA G (0, 1]. 
4 -^fc+1/2 ^k 

This yields the result. □ 

Since Gx^ L\r\ L\, Theorem A. |l] shows that Q^^ {Gx, Gx) G L'^- The equation 

-A^AGA = Qe,(GA,GA) (3.10) 

implies that AGa G L'^- Thus, a bootstrap argument shows the smoothness of Ga. This 
reasoning will not help to find A-uniform Sobolev estimates since the diffusive heating in 
(p.lOh will vanish in the formal limit A — > 0. 
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Theorem 3.6. Assume that e(-) belongs to the class Km for some integer m ^ 2. Then, for 
any k ^ and integer i e [0, m — 1] 

sup IIGaIIh* < oo. 
Ae(o,i] 

In particular, if e{-) belongs to the class with m > 3 one has sup;^g(Q i] ||Ga||l^ < oo for 
any k ^ 0. 



Proof Use induction over \i\ G N. Proposition |3.5| shows that the result is true if \i\ = 0. 
Let then \£\ := s + 1 > be fixed and assume that for any k ^ there exists > such 
that 



max sup II^^^Ga 11/^2 ^ Ck- (3.11) 



Observe that differentiating ^-times Equation ( |1.16p yields 



Multiplying this equation by d^Gx{v){v)'^^ and integrating over we get, as in Proposition 

M 

X^WVd'GxWl^^ [ d^Qe,{Gx,Gx){v)d^Gx{v){vf''dv + {2e + k)\^\\d'Gx\\l2. (3.12) 

k J^3 k 

Fix A; ^ i. One has 

/ d'Ql{Gx^Gx){v)d'G^{v){v)'"^dv^\\d'QX{Gx,Gx)\y WGx\\l'^ ^ 

^\\QUGx^Gx)\\^s+r\\d'Gx\\L2 

fc-i k+^ 

since |£| = s + 1. One estimates the Sobolev norm of Q+^(Ga, Ga) thanks to Theorem ^ 
applied to rj = k — ^. Precisely, for any e > 0, 

l|Q+(GA,GA)||H^..^^G(e)||GA||H|,^,,3l|GA||,.^.^^^3+.||GA||^ IIGaIIh^^, 

+ e(\\Gx\\L^ Jd^GxWL. ,+I|G'a||li JI^^GaH^^ \. 

Using the uniform bounds in given by ( |3.11p together with Proposition and the 
uniform L\ bounds, one notes that there exist a/c, /3fc > such that 

||Q+(Ga,Ga)||h=+i ^a^^ehWd'GxWL^ , VAg(0,1]. 

Therefore, 

I a^Q+(GA,GA)(t;)a^GA(t;)(?;)''=d?;^a,.||a^GA||i2 +e/3fc||9^GA||i2 . (3.13) 
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Regarding the loss part of the collision operator, first note that 



u=0 ^ ^ 



For any / \£\, integration by parts yields 



Q-(d-'Gx,d'-''Gx){v) =\d''Gx{v) 



d'-''Gx{v.)\v-v^\dv, 



^ \d''Gx{v)\ Wd'-'^-'Gxhi 

where, £ — 1 = (^i — 1, ^2, ^3) for any multi-index £ = (^1, ^2, ^3)- Using again the control of 
norms by weighted L^-norms, see inequality ( |2.15p , we get 



Q-(d''Gx,d'-''Gx){v) 



^^T\d''Gx(v)\ \\d^-^-^G 



for some universal constant r > independent of A. From the induction hypothesis (BTTTp, 
this last quantity is uniformly bounded and using Cauchy-Schwarz inequality we obtain 

Q-^{d''Gx,d'-''Gx){v)d'Gx{v){vfUv 

3 

E (^)l|S''GAlli.||fl'0A(t>)llij<C„||8'0Alltj VA€(0,1] 



\u\<\e\ ^ 



\u\<\e 



for some positive constant Gk/ independent of A. Second, whenever u = £ we have ac- 
cording to Proposition p.2| the lower bound 

/ Q-{d'Gx,Gx){v)d'Gx{v){v)^''dv^co\\d^Gx\\l2 . 

Thus, summarizing, inequality ( |3.12p reads 



X^WVd'GxWl^ ^GkAd'Gxhl+akWd'Gxh^ , + e/?^ 
+ Ck\\d^Gx\\Li-co\\d^Gx\\l2 



k+l 



+ {2k^ + k)\^\\d^'Gx\\l2 VAg(0,1]. 



Choose e > such that e/3fc = ^. We note that, after neglecting the gradient term in 
the above left side and bounding all L\ norms by ^^^i norms, there exists some positive 

constant > such that 



^Wd'GxWi^ ^AkWd'GxWL^ +{2e + k)X<\\d'Gx\\i2 VAe(0,1]. 

^ k+l k+l k 



Finally following the proof of Proposition |3.5| , we get that supj^g/g il 11^ ^'aH^.z < 00 for 
any /c ^ i. □ 
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Remark 3.7. In the constant restitution case, uniform regularity estimates where obtained 
using the propagation of regularity and damping with time of singularities for solution to the 
time-dependent problem. More precisely, using the fact that the solution to 

dtf{t,v) = QMJ){t,v) + {l-a)A,f{t,v) 

can be written as f{t, v) = fsit, v) + fnit, v) where fs is smooth and the reminder f^ is small 
in some appropriate norm, see p2[]. Our approach applies to such case yielding a much more 
direct proof of these estimates. 



The proof of Theorem can be easily modified to get an estimate of the difference of 
solutions to ( |1.15p . 



Proposition 3.8. Assume that e(-) belongs to Em for some m ^ 3. For any integer £ G 
[0, m — 1], there exist some positive constant Q > such that 



II^a-GaIIh^ ^QII^^a-GaILi VAg(0,1] (3.14) 
for any G\,Fx e &x. As a consequence, there exists a positive constant Ca > such that 

\\Fx-Gx\\^i,i^^^)^Ca\\Fx-Gx\\L^im.) VaG[0,|], A G (0, 1] (3.15) 
where the weight rria := ma{v) = exp{a\v\). 

Proof We follow the argument of [ |20| , Proposition 2.7] that uses induction and only give 
the details for the initial step £ = 0. Set Hx = Fx — Gx, we aim therefore to control ||i/A||L2 
by Ili/Allii- Notice that Hx satisfies 

Qe,{Hx,Fx) + QeAGx,Hx) = -y'AHx, VA G [0,1]. 
Multiplying this identity by Hx and integrating over yields 

A^||Vi/A|li2 + / Q~iHx,Fx)Hxdv 

= / {Qt,{Hx,Fx) + Qt,{Gx,Hx))Hx{v)dv- f Q-{Gx,Hx)Hxdv. 
From Proposition |3.2| one has 

/ Q-{Hx,Fx)Hxdv;^co\\Hx\\l2. 

In addition, 

/ Q~{Gx,Hx)Hxdv^ f Gx{v)\Hx{v)\dv f \Hx{v,)\\v - v,\dv. ^ G\\Hx\\li 

JR3 Jm.3 J^3 1 

where the constant G depends only on the L°° norm of Gx which is uniformly bounded. In 
order to control the gain operator, split the angular kernel b{s) = into h{s) = hi{s)+b2{s) 
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with 61 (s) := ^l(_i+5^i_5)(s) for some (5 > to be determined latter on. Using Young's 
inequality, see Theorem A. |2| 

/ {Qt^{Hx,Fx) + Ql{Gx,Hx)) Hx{v)dv 

^ C(6i) (IIFaIIl^ + IIGaIIl^) II^Alli; + ^(62) (||Fa|Ili + IIGaII^i) ||^A|li2 . 



Fix e > 0, from the explicit expression of both C{bi) and (7(62) provided by Theorem A. ^ 
one notes that it is possible to choose 6 > such that (7(62) ^ e (recall that b is bounded). 
Summarizing, for any e > 



A^||V/7A|li2 + co||/7A|li2 ^ell^Allia +C{e)\\Hx\\li 



% ^^ii^Aiii. 

where C(e) is a positive constant independent of A. Choosing e = ^ we deduce that 



^^||/7A||i. VAG(0,1] 

-^1 Co 1 



which gives the result for ^ = 0. To extend these estimates to higher order derivatives, one 
proceeds by induction using Theorem p .6| yielding ( |3.14p . To deduce now estimate ( |3.15p 



recall the interpolation inequality given in pO| , Appendix B] (note a misprint in the op. cit. 
where the exponents 1/8 have been replaced by 1/4): For any a ^ and £ ^ 0, there exist 
C{a, i) > such that 

where £q := 8£+^ and b = 12a. According to Proposition |3.4| there exists some c > such 
that ||i?A||Li(m6) ^ clli^AlLj for any A G (0, 1]. Moreover, ( p.l4( ) implies that the H^"-norm 
of Hx can be controlled from above by ||i/A||Li • Combining these estimates we get 

II^a||wm(„„) ^ Cia,£)\\Hx\\];(\\Hxf/,\^^) VA G (0,1] 
which yields the desired conclusion. □ 



3.3. Continuity properties of Q+ as A — ^ 0. We investigate in this section the continuity 
of the gain part Qt{f,g) with respect to the restitution coefficient. We shall prove that, for 
sufficiently smooth functions / and g, the collision operator (/, g) converges strongly 
towards Qf{f, (7) as A in a suitable norm to be specified. 
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Proposition 3.9. For any k ^ 0, there exist some explicit constants C(7, k) and (7(7, k) such 
that 



II (Q+ (/,<?) - Qtif,9)) {v)'\\n-^ ^ C{j,k,a)X^\\fhi^^^^^ \\g Wq^^^^^ (3-16) 
II (Q+ (/,<?) - Qtif,9)) {v)'\\n-^ ^ C{j,k)X^\\f\\Li^^^^ Wshu,^, (3.17) 

Proof. Fix A > 0, a test function (/> G and define il;{v) = {v)'^(f){v). Use tlie weak form of 

QUf,9)-Qt{f^9) to get 

(Qe+(/,5)-Q|(/,5)) {vmv)dv 

= ^ I \u-n\f{v)g{v,)(^{v'^^^) + ^{v?''^)-^{v')-ij{v'S]dvdv,dn 

wliere {v'^^\v^'') denotes tlie post-collisional velocities associated to the restitution coef- 
ficient e\ while {v' ,v'^) denotes the post-collisional velocities for elastic interactions, that 
is, 

v' = V — {u ■ n)n , v'^ = Vi, + {u ■ n)n 
= V — f]x {u ■ n)n , vi^^ = v^, + f3\ {u ■ n)n 
with /3a = Px{\u • n|) = 1±£aM1. Set 

\u ■ n\f{v)g{Vi,) (^tp{v^^^) — ip{v')) dv df^ dn 

and define ^ 

C = C(^i ^) = ^^"^^ — = — (ti • ri)n. 



< 00, 



According to Assumption 

1 — e(r) 
r>0 '^'^ 
thus, for any u, n, A, 

|CK^A^k-nr+^ 

Moreover, for any fixed v,Vi,,n, 

'4}{v^^^)-iIj{v')= ! Vil){v' + sQ- Qds. 
Jo 

These two observations lead to 

Ix ^ '^[^^ A^ f dfdT;j,dn f \u ■ n\'^~^'^ f{v)g{vi,) \\/^l^{v' + s ()\ ds. 

2 Jr3xm3xs2 Jo 

At this point it is important to recognize that for any fixed s G (0, 1] the integral 



/. 

JR-- 



\u ■ n\'~'^'^ f{v)g{v^,) \ V^l^{v' + s()\ dv dvj, dn 



26 R. J. ALONSO & B. LODS 

is just the weak form of the gain part of some peculiar Boltzmann-like operator Indeed, 

set ipiy) = \V'4}{v)\ and V!. = v' + sC, (notice that V!. depends on u,n, A and s) and observe 
that 

Vg = V — j3s {u ■ n)n 

for some parameter 

= PsiW ■ n\) = {l-s) + s^x{\u • n|) G (1/2, 1]. 

Therefore, V!. is in fact a new post-coUisional velocity associated to the above (is- We 
compute for any s G (0, 1], 

u ■ nY''^'^ f{v)g{Vi,) I V'0(w' + s C) | dv^, dn 

\u-n\"<^'^f{v)g{Vi,)ip{V!,)dvdVi,<ln= / ~ {f , g){v)(p{v) dv 

where the collision kernel Bq is given by Bq{u, n) = \u-n\'^^'^ and the restitution coefficient 
eg is such that 13s = Since 

^{v) ^ {v)''\V(p{v)\+k{v)''-^\(t>{v)\ 

^max(l,A:)(|V0(^;)| + |(/)(^;)|)(^;)^ 

one has for any s G (0, 1) 

\u ■ ny^"^ f{v)g{vi,) \ Vil){v' + s C)| duj, dn 



^max(l,A;) / Q% ~ {f , g){v) {\V (t>{v)\ + \(t^{v)\) {v)' dv . 
Jm.3 

As a consequence, thanks to Cauchy-Schwarz inequality 

h ^ V2m8.x{l,k)e^{e)X^\\(t>\\n^ \\QBo,eSf^ a^qds. (3.18) 

It remains to estimate the norm \\Q%g g {f,g)\\i2^ for any s G (0, 1). This is simply done 
using Theorem A. |l| 

In Theorem A. ^ is shown that C(es) only depends on the value at zero of the restitution 
coefficient. Since es(0) = 1 for any s one gets that C{es) is independent of the variable s. 
Thus, estimate ( p.l6( ) follow from ( p.l8[ ). Exchanging the role of / and g in Theorem A. ^ 
gives the second estimate. □ 

We use the equivalence of norms (that follows using Fourier transform) 

\\'^^fn-^ + \Mh-^ = \Mh (3.19) 



valid for any ip € to make Proposition 3.9 stronger. 
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Proposition 3.10. For any i £N and /c ^ there exists C(7, k, i) such that 

\\QtAf,9)-Qtif,9)\\ui 

^C(7,A;,^) (ll/II^M Iblln^+i + II/IIh^+i llffllw^'^ ) 

\ J^fc+7+2 J 

holds for any A G [0, 1]. 
Proof. Since 

VQ+ = Qe+ (V/,5) + Qe+ (/, V5) (3.20) 
and the same is true for Qi{f,g), it suffices to apply Proposition and identity ( ]3.19p 
conveniently to each term to get the conclusion for 1 = 1. Use induction to obtain the 
result for higher derivatives £ > 1. □ 

It is actually possible to extend these estimates to the smaller space L^{ma) with expo- 
nential weights 

ma{v) := exp {a\v\) , u S M^ a ^ 0. (3.21) 
We work with exponent 1 for simplicity even if, as suggested also by Proposition |3.3| , it 
is likely that our results are still valid for general weights of the form exp(a|u|P) with 
^ p < |. The advantage of the following result with respect to the previous one is that it 
involves the derivative of only one of the functions / or g. Precisely, one has the following 
that extends [|19|, Proposition 3.2]. 

Theorem 3.11. There exists an explicit constant Aq S (0, 1) such that for any a ^ there 
exists C(7, a) > 0/or which there holds 

\\Qt,{f,9) - Qtif,9)\\mm^) ^ C7(7,a)A^||/||ii(^^) Mwl'^m.) ^ (0, Aq) 

(3.22) 

and 

\\Qt,if,9) - Qtif,9)\\LH^m.) ^ C'(7,a)A8^||5ki(^„) Il/Ilw^^(,n.) ^A G (0,Ao). 

Proof The proof follows the argument of the analogue [ [T9t Proposition 3.2] where the 
crucial estimate is provided by Proposition A.^ (see Appendix A). Precisely, as in the op. 
cit, for any given v, G M'^, w = v + ^ and a G S^, we define the angle x £ [0, f ], 
by cos X := |o" • t^l- Let (5 G (0, 1) and > 1 be fixed and let bs G W^'°°(-l, 1) such that 
bs{s) = bs{—s) for any s G (0, 1) and 



bsis) 

with moreover 



1 if s G (-1 + 2^, 1 -2J) 
if s ^ {-1 + 5,1- 5) 



^ bs{s) ^ 1 and \bs{s)\ ^ | Vs G (-1, 1). 

Let us define also QR{r) = Q{r/R) with Q{x) = 1 on [0, 1], e(x) = 1 - x for x G [1,2] 
and G(x) = on [2,oo). We define the sets A{5) := {a G S^; sin^ x > 5}, B{5) := {a G 
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; u ■ a ^ (— 1 + 25, 1 — 2(5) or sin^ x ^ f^}- With these notations, for any restitution 
coefficient e(-), we split Qf into 

0+ = 0+ + 0+ + Qt 
where the collision kernels i?i(n, n • ct), i = 0, 1, 2, are defined by 



and 



B2{u,u-a) = bsiu- a)@R{u) — , Bi{u,u ■ a) := — 1a(5) (1 - Qr{\u\)) 



u\ \u\ 
Bo{u, 2 . a) = U (1 _ 6,(2 . a)) e^dnl) + ^ (1 - QrM)) Ia^s)- 

We shall of course apply this splitting to the restitution coefficients e\ and the elastic one 
e = 1 which corresponds to eo- The proof is divided into three steps. 



• Step 1. Estimate for We can prove exactly as in [ ]T9i Proposition 3.2] (precisely, 

using Theorem A. |l[) that for any A G [0, 1) and any 5 £ (0, 1) there holds 



• Step 2. Estimate for We can check without difficulty that \ \T9[ Lemma 3.3] still 
holds true for non-constant restitution coefficient e\ with exponent k in Eq. (3.8)] 
given by /c = (1 — 5/160) 2 independent of A. In particular, reproducing the proof of the op. 
cit. we get that there exists a constant C > such that: 

l|25„e.a5)|Li(„„)^^ll/llL;(-»)NlL}(™„) VAG[0,1],<^G(0,1), R>1. (3.24) 

• Step 3. Estimate for the difference — ^: The crucial point is now to estimate 
II Qfia.e^ 5') - Qb2,iU\ 5) IIli (m,) ^nd, as already mentioned, we shall resort to Proposition 
A. ^ given in Appendix A. Precisely, let cj) £ L"^ and '4){v) = ma{v)(l){v). Notice that the 
collision kernel i?2(^i, S-o") satisfies the assumption of Proposition A. |2] since Supp^^ C (—1+ 
(5, 1 — 5). Applying this Proposition to the restitution coefficient e\ (with fixed A G (0, 1]) 
one sees that there exists Ce;^ > such that 

1 



+ 2T+%(eA) ds Q+ {f,h)\^P{v)\dv 

Jo Jr3 ^T'^s 

where h{v) = g{v) + |V5(?;)| while the kernels and are given by 

B^{u, u ■ a) = B2{u, u ■ a)\u\'^ , B^{u, u ■ a) = max{B2{u, u ■ a), \VuB2{u, u ■ a)\)\u\'^~^'^ 

and, for any s G [0, 1], e^(-) is a given restitution coefficient with in particular e^(0) = 1 
for any s, A. One estimates these two integrals using Theorem A. [l[ Precisely, by Holder's 
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inequality 



QtAf,9)mv)\dv^\\Q+,{f,g)ma\\Li 



while, for any s G (0, 1) 



Now, one notices that 

ina{vs) ^ ma{v) ma{v^) 



and 



rriaiv'i) ^ ma{v) ma{Vi,) 



where v'^ and denote the post-collision velocities associated to the restitution coefficient 
and e = 1 respectively so that \\QQ_^^i{f,g)ma\\L^ ^ WQ^^^ii'^a f,rnag)\\Li and 

Since 0R(|n|) = whenever > 2R, one has eij(|n|)|'u|'^ ^ {2R)'^ for any n G and 
there exists an universal constant ci > such that 



(3.25) 



To estimate ||Q± _^ (ttiq /, ma h) one only notices that the kernel u ■ a) can be 
estimated by 



B.y{u, u ■ a) ^ 







[a 


+ 1) i„i + 1_ 







bs{u ■ a)ER{\u\) \u 



7+2 



for some positive mapping such that = whenever \u\ > 2R; the factor 

coming from the derivative of while the term 6"^ comes from that of 65. Then, using 
as above Theorem A. |^ and because e^(0) = 1 is independent of s £ (0, 1), one gets the 
existence of an universal constant C2 > such that 

"^i,e>' {rnaf,mah) 11^1 ^ C2 (^-^ + WrUahW^i 

or, equivalently, 

gA ("T-a/^Wla/i) IIli ^ C2 ^"^ + ^^^^ ||"la / HlJ ||"laC?||^l,l ■ (3.26) 

Finally using the fact that i^{e\) ^ X"/ i^[e) while, as noticed in Remark A. |], Ce^ ^ C3A''' 
for any A G (0, Aq] for some constructive Aq > and some positive constant C3 > 0, we 
finally obtain, combining ( |3.25p and ( |3.26| ) that 



^ CA^ ( i + i ) 



L\(m,)\\9\\w\'\^a) VAg(0,Ao] (3.27) 
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for some positive constant C > 0. Collecting estimates ( |3.23| )-( p.24t )-( p^27[ ), we finally 
get that there is some positive C > such that 



\Qt{f^9)-Q-^{f,9)\\LHm.) 



^ C { 6 + + 



VA G (0,Ao), (5 > 0, R>1. (3.28) 



Then, choosing 5 and R> 1 such that 





for some p > 0, one sees that necessarily p = and R^X^ = X^p. This gives the 

conclusion. One proves the second estimate exactly in the same way. □ 

Notice that, increasing the pol5Tiomial weights in the various norms of / and g, we can 
get an optimal control rate A'*^. 

Corollary 3.12. There exists some explicit Aq € (0, 1) such that for any a ^ there exists 
some explicit constant C{-f, a) > 0/or which there holds 

\\Qt,if,9) - Qtif^MLHm.) ^ ^^(7,a)A^||/ILi(™,) \\g\\wi'\^.) e (0,Ao) (3.29) 



and 



VAg (0,Ao) 



where k = j + 



Proof. The proof follows the lines given for Theorem [3.11 . Let us explain the small changes. 
Bounding directly 6r by 1 allows to replace Estimate ( |3.25p by 



In the same way, for some given l</c<7 + 3tobe determined later, one can replace 
Estimate ( |3^ by 



With such a choice, (3.27) becomes 



+ C^A^II/llL^,,(™.)lbllL^^,(™„)VAG(0,Ao] 
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and, collecting all the estimates as above we get 

\\QUf,9)-Q''{f,9)\\LHm.)^ 

(P7+3— fc\7 \ 

VA G (0,Ao), 5>0, R>1 

where s = max{k, 7 + 1). One looks now for A; G (1,7 + 2) for which it is possible to choose 
6, R> 1 such that 



and we get that, necessarily, /c = 7 + 4^. In this case, R^^'^^^X' = A^'^ and we obtain finally 



4. Uniqueness 



We are now in position to prove the uniqueness of the solution to ( |1.15| ) for sufficiently 
small A, that is, there exists A^" > such that for any A G (0, A^) the stationary problem 
( |1.15 ) admits a unique solution Gx with unit mass and vanishing momentum. The strategy 



of proof has been sketched in the Introduction and we shall refer to Section |1.3| for the 
main steps of the proof. In particular, a crucial point consists in proving and quantifying 
the convergence of {G\)\ towards an universal limit M. This is the object of the following 
paragraph. 

4.1. The limit A — ;> 0: non quantitative version. Using the continuity properties, specifi- 
cally Theorem [3.11| , and a compactness argument, we establish a first convergence result, 
non quantitative in the sense that no rate of convergence is provided. 

Theorem 4.1. Assume that e(-) belongs to the class Em with m > |. For any k ^ and 

i G [0,m — 1], one has 



lim ||Ga - A^IIh^ = 

A-J-O 



where M is the Maxwellian 



3 / w'^ 



M{v) = (2^G)-2exp( 



The Maxwellian's temperature Q is given by 



2 



sJ^^Y' C4.1) 

\ a 2 2 m3+-y J 



where m^^^ is the (3 + ^)-th moment of a normalized Gaussian 

.if ( kn, ,3+. , ^3±.r(3 + i) 
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Proof. The proof is divided in several steps and essentially based upon a compactness ar- 



gument through Theorem 3.6 



• First step: compactness argument. Let us choose m — 1 ^ £ > | and A;o ^ 1 in 
the Theorem |3.6| . It clearly exists a sequence (A„)n with A„ — ;> and Go G such 
that (Ga„)„ converges weakly, in H^^, to Go (notice that, a priori, the limit function Go 
depends on the choice of £ and ko). Using the decay of Ga guaranteed by the polynomially 
weighted Sobolev estimates, we can prove thanks to a simple localization argument (and 
using compact embedding for Sobolev spaces) that the convergence is actually strong in 

for any k < ko 

lim 1 1 Ga„ - Go 1 1 Hi =0. (4.2) 

n— >oo k 

Indeed, since supAg(o,i) ||Ga — GoUh'^ < any fixed ^ k < ko and any e > 0, there 

is > large enough such that 

sup ||Ga - Go||h1(bs,) ^ ^ (4.3) 

Ae(o,i) 

where Bji = {v G , \v\ ^ R} and Bf^ its complementary. Let Ga„ and Go denote the 
restrictions of Ga„ and Go to the ball Bfi. Since ^ > |, according to Rellich-Kondrachov 
compactness theorem [|l]. Theorem 6.2, p. 144], the embedding il^{Bji) ^ il^{Bii) is 
compact so that there is a subsequence of (GA„)n that converges strongly to Go in H^(i?ij). 
Since Go is the unique limit of all subsequences, it is actually the whole sequence (Ga„ )n 
that converges to Go in H^(i?R). Combining this with ( ]4.3p yields (4.2). 

• Second step: identification of the limit Gq. Let us prove now that the above limit Go is 
actually a Maxwellian distribution with temperature 6. To do so, one uses ( |1.15p to get 

||Qe,(GA,GA)||L2 = A^||A,Ga||l2 VA > 



and, since supAg(o,i] ||A„Ga||l2 ^ supAg(o,i] ||Ga||h2 =: Go < oo according to Theorem [3T6 
we get 

||Qe,(GA,GA)||L2 ^GoA^ VAe(0,1). (4.4) 
Now, from the identity Q^iGx,Gx) = Qe^{Gx,Gx), one has ||QeJGA,GA) - 
Qi (Ga, Ga) ||l2 = II Q+ (Ga, Ga) - Qf (Ga, Ga) 11^2 so that, 

l|Ql(GA,GA)||i2 ^ ||Q+(GA,GA)-Q^(GA,GA)||L2 + ||Qe,(GA,GA)||i2. 

Combining the above estimate ( |4.4p with Proposition |3.1C1| we get 

||Qi(Ga,Ga)||l2 ^GoA^ + GiA^||Ga||^i 



for some positive constant Gi > independent of A. Using again Theorem |3.6| , we get that 
there exists some explicit constant G2 > such that 

I|Qi(Ga,Ga)||l2 ^G2A^ VAg(0,1]. 

In particular the sequence (GA,Jn constructed in the first step satisfies 

lim ||Qi(GA„,GAj||i2 = 0. 
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Since Ga„ — ^ Go strongly in Lf, we get easily that 

Qi(Go,Go) = 

i.e. Go is a Maxwellian distribution. By conservation of mass and momentum, we get that 
Go has unit mass and zero momentum and it remains only to determine its temperature 
9. To do so, we shall use equation ( |1.16| ) and Lemma A. ^ With the notations of Lemma 
A. ^ equation ( |1.16p writes X\{G\) = 6 for any A G (0, 1]. Applying Lemma A. |^ with 
fi = 91 = Ga„, f2 = 92 = Go (with for simplicity 6 = 1) and estimating the weighted 
-norms by L^-norms (using equation ( |2.15p with = ^ for instance) we get that, for any 
e > there is no ^ 1 such that 

|Xo(Go)-6KGi||GA„-Go||i2 +G2e Vn ^ no 

5+7 

for some positive constants Gi, G2 > independent of n where we used that A;o ^ 6 + 7 

and the uniform estimates on ||Ga|| 7-2 . Letting n go to infinity, we get that Xo(Go) = 6. 

5+7 

Therefore, 

6 = C^ [ Go{v)Go{v^)\v-v^f+^dvdv^ = C^e^ [ M{v)M{v^)\v-v^f+^ dv dv^ 

where M(-) is the normalized Maxwellian M{v) = 7r~i exp Some algebra yields 

6 = G-y0^22m3_|_^ from which we deduce (4.1), and thus. Go = M. 

• Final step: convergence of the whole net (Ga)a- We conclude the proof by showing that 

lim ||Ga - A^IIhi = 0. 

A->-0 

Argue by contradiction assuming this does not hold. Then, there exists eg > and a 
sequence (A„)„ converging to zero such that 

\\Gx,. - M\\ni ^ eo VnGN. 
We just proved above that (Ga„ )n admits a subsequence (Ga„ )j converging strongly in 
to Ai. Therefore, 

eo ^ ||Ga„, -A^IIhi 0, 

which is a contradiction. This proves that the full net (Ga)a converges to M strongly in any 
H^. We proceed along the same path (using a version of Rellich-Kondrachov Theorem for 
higher-order Sobolev spaces) to prove that the convergence actually holds in any weighted 
Sobolev space H^, k ^ and ^ G [0, m — 1]. □ 

This convergence in Sobolev spaces can be extended easily to weighted L^-spaces with 
exponential weights. Recall that, for any a ^ 0, we denote 

ma{v) = exp(a|t'|), v G M^. 

Corollary 4.2. Assume that e(-) belongs to the class Em with m > |. For any a ^ and any 
k>Oit holds 



lim||GA-A^|Li(^^) = 0. 
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Proof. Taking ^ > | in the above Theorem, observe that by classical Sobolev embedding 

lim IIGa -M\\roc = 0. 

The proof follows then using interpolation. First, observe that the convergence holds in 
exponential weighted L^-spaces 

\Gx{v) - M{v)\'^mb{v)dv !^Cb\\Gx- M\\l-- ^0 as A^O 



where Cb := supAg(o,i] ||Ga — -A^|lLi(m5) finite for any 6^0 thanks to Proposition 
Then, using Holder's inequality, for any b,a > 

The last integral in the right side is finite provided b > 2a, therefore, the convergence 
follows from the convergence just proved. □ 

4.2. Uniqueness result. On the basis of the above convergence result, we are in position 
to apply our general strategy as explained in Section 13. Recall that, for any given A G 
(0, 1] and Gx, Fx G ©a, we set 

Hx = F\ — Gx- 

We have to determine Banach spaces X and y for which the estimates ( ]1.20| ) - ( |1.24| ) hold 
true. The analysis of the previous section suggests the choice 

X = L\ma), y = L\{ma) 



for some exponential weight niaiy) = exp(a|ti|), a ^ 0. Indeed, Proposition p.ll and 
Corollary ^ already ensure that (|L2lt) and ( [lT23l ) are fulfilled. The fact that qi.20p 
stands is a classical property of Boltzmann operator with hard-spheres interaction (see 
Theorem 12]). Since 

||Ai?A|U ^ II^^A||w2'i(ma) 

one sees that estimate ( I.22p holds because of Proposition (precisely, inequality 
( |3.15[ )). Now, property ( |1.24 ) is a consequence of the following spectral property of ^i, 
first established in [21]. 

Proposition 4.3. The spectrum of the linearized operator in X (with domain ^(Jfi) = y) 
has the following structure: 

(1) is a simple eigenvalue of associated to the null set 

J^i^i) = Sp&n{M,viM,V2M,V3M,\v\^M); 

(2) the continuous spectrum of is given by (— oo, —uq] where 

vq = inf / \v - Vi,\M{Vi,) du^; 



(3) the non zero eigenvalues of J£\ are all negative and can accumulate only at —uq. 
Consequently, ^\ admits a positive spectral gap u > 0. 
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In particular, if 

x = {fGX;[ fdv= [ vfdv= [ \v\^f{v)dv = o}, y = ynx 

JrS Jr3 J]fj3 

then ,yV{J^i) r\ y = {0} and Jfi is invertible from y to X with explicit estimates for 



-^1 ^Ili^^S? ■ Co^^^Qu^^tty! inequality ( 1.24 ) holds true with cq = ^ 



Remark 4.4. The proof of the above proposition can be seen as a consequence of some gen- 
eral comparison principle that asserts that the linearized collision operator enjoys the same 
spectral properties in X and in the largest Hilbert space T-L = L'^{M.^^). A simple proof of 



Proposition4.3 can also be recovered from 



The difference H\ = Fx — G\ does not necessary belong to y since we do not know 
a priori that Gx and Fx share the same kinetic energy. Consequently, we need a slight 



modification of the strategy developed in Section |1.3| to state our main result, regarding 
uniqueness of the steady state. 

Theorem 4.5. Let e(-) belong to the class Em for some integer m ^ 4. There exists G (0, 1] 
such that 



6x = \Gx € lI;Gx solution to dTTsl ) with / Gx{v)dv = land / 



vGxiv)dv = 



reduces to a singleton for any AG [0, A^). 



Proof We explained in the previous paragraph that the estimates ( ]1.2Clp , ( |1.21p , ( |1.22| ) 



( |1.23D and ( |1.24D of the general strategy are fulfilled with X = L^(m,a) and y = Ll{ma) 



for any a ^ 0. Let us fix e > and reproduce the computations of Section ^]3|. It follows 
that there exists Aq G (0, 1) such that 

\\^i{Hx)\U ^ e\\Hx\\y VAg(0,Ao). (4.5) 
Let us now introduce the following lifting of the operator into an invertible operator 

A : h^Ah= {Aih; A2h) G M x ^ (4.6) 
where the second component A2h = while the first component is defined by 

Aih = 2Io{M,h) = 2 M{v)h{v^)Co {\v - v^\^) dvdv^. 

JR3xIR3 

We refer to the Appendix A for notations. Since Gx and Fx share the same mass and 
momentum, one deduces from Lemma A. |^ that 

Hx = Fx-Gx = A-UHx = A-^ {Ai{Hxy,A2{Hx)) 
with an explicit estimate of the norm ||^~^||. Since A^^ maps M. x X toy, we get 

ll/^Ab^ P^'ll m^x(^\Ai{Hx)\; P2(^a)Iu)- (4-7) 
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Still using the notations of Appendix A, one has readily 

+ (lo{Fx + Gx, Fx - Gx) - IxiFx + Gx, Fx - Gx) 
Now, it is clear that 

\UM - Fx,Hx)+Io{M - Gx,Hx)\ ^ C, f\\M - Fxh^ + \\M - Gxhi^ ) 11^^11^^+ 

\ 3+7 3+7/ 3+7 

^ G^ {\\M - FxWx + \\M - GaIU) II^Ab VA G (0, 1]. 

Consequently, applying then Lemma A. |^ with fi = fi = Fx + Gx and gi = §2 = Fx — Gx, 
we get the existence of constant C > and some Ai G (0, 1) such that 

\Ai{Hx)\^C{e+\\M-Fx\\x + \\M-Gx\\x)\\Hx\\y VAg(0,Ai). (4.8) 

In particular, on the basis of ( |1.23p , there exists A2 G (0, Ai) such that 

\Ai{Hx)\<.2Ge\\Hx\\y VAg(0,A2). 

Using now ( f4.5p together with ( 14.7D , we deduce that 

\\Hx\\y ^ emax(l,2C)p-^|| \\Hx\\y VA G {<d,\^) 

where A^" = min(Ao, A2). Taking e > small enough yields therefore the desired unique- 
ness: Fa = for any A G (0, A^). □ 

4.3. Quantitative version of the uniqueness result. We derive is this section a quantita- 
tive version of the Theorem 4J. which shall result in a quantitative estimate of the above 
parameter A^. 

Proposition 4.6. Let e(-) belongs to the class with m ^ 4. Assume moreover that there 
exist two positive constants a, b > and two exponents 7 > 7 > such that 

|e(r) — 1 + ar'^l ^ b r'^ for any r ^ 0. (4.9) 

For any a ^ 0, there exist some explicit Aq G (0, 1), cq, ci > and exponent a = min(7, 7 — 
7) > such that the estimate 

I|Ga-A^||li(™„) ^coA° + ci||GA-A1|lii(^^) VAg (0,Ao) 

holds for any Ga G ©a- 

Proof We apply a slight modification of the proof of Theorem |4.5| where, instead of esti- 
mating the difference of two solutions to ( 1.15p , we estimate the difference Gx —M.. Recall 
that A is the lifting operator given by ( |4.6| ). Thus, 

Gx-M = A-^A{Gx -M)= ^"^ (A(Ga - M);A2iGx - M)) 

where the norm of ||^^^|| is explicit. In particular, since A^^ maps M x A" to 3^, we get 

||GA-A^b ^ \\A-'\\max(^\Ai{Gx-M)\;\\A2{Gx-M)\\;,y (4.10) 
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Let us estimate separately the two terms Ai{Gx — M) and A2{Gx — M). In the one hand, 

Ai{Gx -M) = 2Io{M, Gx-M)= lo{M - Gx, Gx - M)+ 

(^o{Gx, Gx) — Tx{Gx, Gx) 

where we used the fact that Zq{M,M) = IxiGx, Ga) = 6. Now, it is clear that 

\IoiM-Gx,Gx-M)\i^C4Gx-M\\l, . 

3+7 

Moreover, according to Lemma A. |^ and under assumption ( 14.9D , 

^o{Gx,Gx) — lx{Gx,Gx) ^ C'oA"||Ga||?i 

-^3+7+7 

for some explicit constant Co > and exponent a = min(7, 7 — 7) > 0. Therefore, since 
supAg(o,i] ||Ga|||i , with the notations of the previous section 



^3+7 



\Ai{Gx-M)\^CiX'' + C4Gx 
for some positive constant Ci > 0. On the other hand, 

A2iGx -M) = Qi{Gx-M,M-Gx) + QiiGx,Gx) 
= Qi{Gx - M,M - Gx) + {Qi{Gx,Gx 
= Q,{Gx - M,M - Gx) + {Qi{Gx,Gx 



VA G (0, 1] 



(4.11) 



Qe,{Gx,Gx)) + QedGx,Gx) 
Qe,{Gx,Gx))-X^AGx. 



Therefore, using Corollary |3 .121 , there exist an explicit Aq G (0, 1) and constants (72,6*3 > 
such that 



\A2iGx - ^ C^WGx - + C^2A^||Ga||z.i 



7+^ 

+ A'''||GA||w2.l(ma) 



^1,1 

VAe (0,Ao) 



(4.12) 



Using interpolation, similar to the proof of ( l3.15| ), and Proposition |3.4| we obtain 

sup ( \\Gx\\wi,i + \\Gx\\w^,i(m^) ) < 00. 
Ae(o,i] V 1+^ / 

Hence, inequality ( [4.12^ reads 

\\A2{Gx-M)\\mm.) ^ C'iIIGa + C3A^ VA G (0,Ao) 

for explicit constants Ci , C3 > 0. Combining this estimate with ( |4.1ip and ( |4.10P yields 
the desired conclusion. □ 



Theorem 4.7. Assume that e(-) satisfies (4.9) and belongs to the class with m ^ 4. Fix 
the exponential weight ma{v) = exp(a|7;|) with a ^ 0. There exist an explicit X* G (0, 1) and 
constant c > such that 



l|G'A-A^||LHm„)^cA° 
where Gx G ©a cind a = min(7, 7 — 7). 



VA G (0, A* 
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Proof. The proof follows from the Proposition 4^ and the non quantitative convergence 
Theorem |4.1| . Indeed, recall the estimate 



l|GA->l|lLl(,n„) ^CoA" + Ci||GA->I||ii(^^) VA € (0, Aq) (4.13) 

for some explicit constants co,ci > 0. Then, since liniA^o ||Ga — M-W^i^^^^-^ = 0, there is 
some a priori non explicit A* G (0, Aq) such that 

ci\\Gx-M\\li^^^) ^ \ VAg (0,A*). (4.14) 
Therefore, estimate ( |4.13p becomes 



l|GA-^^||Ll(n^J ^2coA" VAg(0,A^). (4.15) 
This gives a posteriori an explicit estimate for A* since the optimal A* will be the one for 



which ( |4.14| ) and ( l4.15p are identity which yields the estimate A* ^ (4coCi) q. Since all 



the parameters cq, ci, a happen to be explicitly computable we get the result. □ 

Remark 4.8. We wish to emphasize here several points about our approach. First, recall that 
in the case of constant restitution coefficient, the approach of pOj ] yields directly quantitative 
results. This M/as possible thanks to a clever application of the Cercignani's conjecture for the 
elastic Boltzmann operator derived in p6|]. This allowed to compare the entropy dissipation 
functional and the distance to a given Maxwellian distribution, more specifically, the distance 
to the elastic limit. The disadvantage of this approach is that requires pointwise exponential 
lower bounds and high regularity for the associated steady solution. Such lower bounds are 
related to the spreading property of the collision operator and their technical extension to the 
case of non-constant restitution coefficient is not trivial. In contrast, the strategy here does not 
uses entropy techniques at all, thus, it does not require neither pointwise lower bounds nor 
regularity assumptions. This makes it well-suited for problems in which no regularity of the 
steady solution is available, see 



It is easy to deduce explicit estimates for the parameter A^^ in Theorem |4.5| under the 
above assumption ( |4.9p on e(-): 



Theorem 4.9. If e(-) belong to the class Em for some integer m ^ 4 and satisfies ( K.9D , there 
is an explicit parameter A^ G (0, 1] such that 6\ reduces to a singleton for any A G [0, A^^). 



Proof Recall that the only non quantitative part in the strategy described in Section 1.3 
was the convergence rate of Gx towards M. It is made explicit now thanks to Theorem 4.7 



and, resuming the above strategy one gets that there exists some explicit Co > such that 

\\^iiHx)\\x < CoX'^WHxWy VA G (0,A^) 
where Hx = Fx — Gx with Fx,Gx G ©a and A* is the parameter in Theorem |4.7|. Recall 



that A* can be estimated from below in an explicit way. Using Theorem 4.7 one can replace 
( |4.8| ) in the proof of Theorem by the following quantitative estimate 

\Ai{Hx)\^CiX''\\Hx\\y VAg(0,A'^) 
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where Ci > is an explicit constant. Then, resuming the proof of Theorem yields 



A < At = minfA*,Co"^/"| . □ 



||i^Ab^C2A"||i/Ab VAg(0,A*) 
where C2 = max(Co, Ci)||^^^|| > is explicit. We see therefore that H\ = provided 

The above uniqueness result in the quasi-elastic limit A — )• translates to a weak ther- 
malization uniqueness result. 

Theorem 4.10. For any restitution coefficient e(-) belonging to the class with m ^ 4 and 
satisfying ( |4.9p , there exists some explicit > such that for any fi ^ fi\ there exists a 
unique solution F to 

Qe{F,F) + fiAF = 
with /jg3 F{v) dt; = 1 and /j^s vF{v) dv = 0. 



asserts 



Proof The proof is a simple consequence of our scaling choice. Indeed, Theorem 4.9 
that, for A < At the steady problem 

admits an unique solution with unit mass and vanishing momentum. Performing the back- 
ward scaling F{v) = X^^Gx {X^^v) one gets that there exists an unique solution F with 
unit mass and vanishing momentum to the problem 

Qe{F, F) + X^+'f AF = 
whenever A < At. This clearly yields the conclusion with /xt = (At)^"*"^ . □ 



Appendix A: Properties of the collision operator 



We collect in this Appendix some facts about the Boltzmann collision operator important 
in their own right. Some of the properties of Qg that we will establish here are known 
and some others new. We shall consider a collision operator with more general collision 
kernel than the hard-spheres case considered in the paper, more precisely, a collision kernel 
B{u, a) of the form 

B{u,a) = <i>{\u\)b{u-a). (A.1) 

The kinetic potential <!>(•) is a suitable nonnegative function in and the angular kernel 
b{-) is assumed in L^(— 1, 1). The associated collision operator QB,e is defined through the 
weak formulation 

/ QB,e(/,/)WWdt; = ^ / f{v)f{v,)AB,em^,^*) <i^*<iv (A.2) 
for any test function ip = iIj{v) where 
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with f v'^ are defined in ( |1.6| ). For any fixed vector u, the angular kernel defines a measure 
on the sphere through the mapping cr e i-> h{u ■ u) e [0, oo] and we will assume it to 
satisfy the renormalized Grad's cut-off assumption 

||6||^i(s.) = 2^||6||^i(_,^,) = l. (A.3) 

For technical reasons, we shall also assume that 

6 : X G [—1, 1] I — > b{x) = b{x) + h{—x) is non decreasing. (A. 4) 

A particularly relevant model is the one of hard-spheres corresponding to <^*(|n|) = \u\ and 
h{u ■ a) = 1 /47r. For this particular model we shall simply denote the collision operator 

QB,e by Qe- 



A.I. Convolution-like estimates for ^. We begin by recalling some of the regularity 
and integrability properties of the gain part established in and [^. We start first 
with Young-like estimates in with > 0. 

Theorem A. 1 (Alonso-Carneiro-Gamba [0]). Assume that the collision kernel B{u,a) = 
^{\u\)b{u ■ a) satisfies ( |A.3p and <!>(•) G L'^f.for some k € M.. In addition, assume that e(-) 
fulfills Assumption [IT1| . Let 1 ^ p, g, r ^ oo with 1/p + 1/q = 1 + 1/r. Then, for any a ^ 0, 
there exists Cp^r,a,k{b) such that 



\\QtM^9)\\Lr^^Cr,p,e.Ab)\\nLZ 

where the constant Cr,p,a,k{b) is given by 

r' 

-3/2r' \ 7 



Cr,p,a,k(b') — C/j , 



-1 



1 



b{s)ds 



1 + s 



+ (1 - /So)- 



.1 



3 

'2r^ 



bis)ds 



for some numerical constant Ck.a,r independent ofb and e(-) and with /3o = /3(0) = li|M 
Theorem A.|l| has been modified in ^] to provide bounds with rj ^0. 



(A.5) 



Theorem A. 2. Assume that the collision kernel B{u, a) = (t{\u\)b{u ■ a) satisfies ( |A.3| ) and 
$(•) G L^^, for some k E In addition, assume that e(-) fulfills Assumption Then, for 
any 1 ^ p ^ oo and 77 G M, there exists Cn^p^B) > such that 



where the constant p ^(5) is given by: 

Cr,,p,k{B) = Ck,n,p l{r],vM 
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with a constant Ck^n,p > depending only on k, r] and p. Furthermore, the dependence on the 
angular kernel is given by 

j{rj,p,b) = J^^^-^j bis)ds, (A.6) 

where 1/p + 1/p' = 1 and 77+ is the positive part ofij. Similarly, there exists C,j p ^.(i?) > 
such that 

25,.a5)|L^c,,,,,(i?)ibiLi_, 11/11^ 

where the constant p ^(i?) is given by 



for some constant Ck,n,p > depending only on k, rj and p. The dependence on the angular 
kernel is given by 

I _£±!l+ 
7(^,P,ft) = /^(^ + (l-/3of^) '^b{s)ds 

where 1/p + 1/p' = 1 and /3o = /3(0) = 

Corollary A. 1. Assume that the collision kernel B{u, a) = ^{\u\)b{u ■ a) satisfies ( |A.3D and 
$(•) G L°°i^ for some k e M.. In addition, assume that e(-) fulfills Assumption Then, for 
any 1 ^ p ^ 00 and rj e M., there exists a numerical constant Ck,n,p > (which does not 
depend on B {■,■)) such that 

' L*J ^ ' " l>7 + '=l + l')l '7 + *: 

A.2. Useful change of variables for non constant restitution coefficient. We establish 
here several changes of variables that are useful for the study of the continuity properties 
given in Section |3.3| . 

Definition A. 1. A restitution coefficient e(-) : r i—?- e{r) e [0,1] is said to belong to the class 
Co if e(-) satisfies the following: 

(1) The mapping r G M+ i-)- e(r) G (0, 1] is absolutely continuous and non-increasing. 

(2) The mapping r G M+ i-)- -deir) ■= r e(r) is strictly increasing. 

(3) e(0) = 1. 

Moreover, for a given 7 > 0, we shall say that e(-) belongs to the class Cj if it belongs to Cq 
and there exists a > such that 

e(r) ~ 1 — ar"' as r ~ 0. 

Remark A. 1. Recall that if e{-) belongs to the class Cj then 

n f \ 1 ~ e(r) 

t^(e) := sup < 00. 

r>0 f'^ 
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Lemma A. 1. Define /3e{r) = and the mapping 

rje : r £ i — > r/3e(r). 

Then, e(-) belongs to the class Cq if and only if rie{-) is strictly increasing and differentiable 
with 

T 1 71 \V] 

- ^ T?e(r) ^ r ; - ^ ^e(^) ^ ~ — f^^ '^"J r > and rj'^{0) = 1. 
Equivalently, the inverse mapping ae(-) ofrie{-) satisfies 

r ^ ae{r) ^ 2r ; ^ a' (r) ^ 2 /or any r > and a' (0) = 1. 

r 

Lemma A. 2. Lemma 2.3] For any cr e and 5 G (0, 2) de/ine the cone 

ns = ^s{(j) = \ueR^\{0} ; u-a>5-l\. (A.7) 



3 , / N + \u\cr 



Define the mapping cls 



Then, is a C°°-diffeomorphismfrom Qs onto Qs* where ^* = 1 + y | with Jacobian 

Ja{u) = l{l + u ■ a) . Its inverse mapping (p^ = <l>^^ is given by (pa{w) = 2w - — — a. 

w ■ a 

With the notations of Lemma A.|l| and Lemma A.|2] we can establish the following change 
of variables formula which generalizes [|17[ Prop. 3.2]. 

Lemma A. 3. For a given restitution coefficient e{-) inthe class Co, one defines the mapping 

Then, for any 6 > 0, He is a C°°-dijfeomorphism from $7^ onto itself with Jacobian Je{\z\) 
given by 

Je{g) = ^{l + KiMQ)))l3!MQ)) ^Q^O. (A.8) 
The inverse mapping vrg = is given by 

T[e{z) = Z 



\A /3e(«e(k|))' 

If one combines the two applications He o we get the change of variables 

Z = Pe{\^a{u)\)^^{u) 

which is a C°°-dijfeomorphismfrom onto Q^*- Its inverse mapping is given by 
with Jacobian given by Jcr{Ce{z))Je{z). One has 

Ce{z) = IIe{z)iPa{z) with fl^iz) = ^^^IMl. 
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Proof. The properties of lie are proven by direct calculations and noticing that ii z = Ile{w) 
then 

z = Q and = ae(|z|). 

With this identity, one can computes the Jacobian of the transformation passing to polar co- 
ordinates. The final expression of Ce(^) is immediate after noticing that ifcr{rw) = r(pcr{w) 
for any r > and any w G M^. □ 

Remark A. 2. Observe that for e(-) belonging to Cq, since '^'^{r) = re'{r) + e(r) ^ 1/or any 
r ^ and f3e{r) G [i, 1], one has the universal bound 

l^Je{Q)^l 0. (A. 9) 

o 

Lemma A. 4. Let e(-) be a restitution coefficient in the class Cq and let s € [0, 1]. Then, there 
exists a restitution coefficient es{-) belonging to Cq such that 

1 + S ifleiz) - 1) = Hesiz) VzeM^ 

where has been defined in Lemma A.^. 

Proof Define ij.{z) = 1 + s{ne{z) — 1) = (1 — s) + s^ie{z) and recall that 

, . ae.{\z\) 



\z\ 



In order to prove that there exists es(-) in the class Co such that fi = fig^, thanks to Lemma 
A. |l], it suffices to prove that the mapping a : r i-)- (1 — s)r + sae(r) satisfies 

r ^ a(r) < 2r ; ^ a'(r) ^ 2 for any r > and a'(0) = 1. 

r 

Since ae satisfies all these properties, it follows that the same is true for a. □ 

The following proposition is reminiscent of the so-called cancellation Lemma for the 
classical Boltzmann operator |25| |. 

Proposition A. 1. Let e(-) be a given restitution coefficient belonging to the class Cq and let 

13{u,u-cj) = e{\u\)b{u-a) 

be a given collision kernel with @{r) ^ and b{s) = b{—s) with Supp b G [—1 + 6,1 — 6] for 
some 5 > 0. Let Qg ^ and Qg ^ denote the positive part of the collision operator associated to 
B with restitution coefficient e(-) and elastic interactions respectively. For any test function ip 
and any given f, g, one has 



f{v)dv / da I iIj{v + z 



tp dv = 

1 



Fv,a iCeiz)) - F^^^i^^iz)) 



Where F„,.(n) = e{\u\)b{u • a)g{v + u) with b{u ■ a) = = u G 



Je{\z\) 

Ja{u) 



dz (A. 10) 
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Proof. We set for simplicity 



Qteif^9){v)Hv)dv and h 



Thus, 



13{u, u ■ a) f {v)g{Vi,)ip{v'^) dv dv^ da 



where u = v — v^, and vL 



^u^-j is the post-collisional velocity associated 



to e(-). In particular, the change of variables v^, ^ u yields 

1 



f{v)dv / da B{u,u ■ a)g{v — u)'il^{v'^) dv duda. 



§2 

The change of variables u — t- —u in the last integral gives 
1 



fiv)dv da e{\u\)b{u-a)g{v + u)'ilj{v + P{\^^{u)\)<^^{u))du 



where we used that, for fixed a, the support of b is included in [—1 + 6,1 — 5] so that the 
variable u belongs to the cone fls defined by ( ]A.7[ ) . With the notations of Lemma aJs], we 
perform the change of variables z = He o ^cr{u) in the previous integral to get 



f{v)dv / da il){v + z) 



Je{\z\) 



F„,, iCeiz)) dz 



(A. 11) 



where Fy^aiu) = 0(|u|)6(n • a)g{v + u). In the same way, for the particular case of elastic 
interactions (i.e. for e = 1) since C,i{z) = fa{z) and Ji(|2;|) = 1 one simply has 



h 



1 



f{v)dv da il;{v + z)Fv,a {^a{z)) dz 



which clearly gives ( |A.10[ ) . 



□ 



Proposition A. 2. Under the assumptions of Proposition A.^ if e(') belongs to the class C^for 
some 7 > then there exists C'e > such that 



^Cef Ql,{f,g)mv)\dv 

+ 2^+6^^(e) [\s [ Q± {f,h)\i;{v)\dv 



Qteif^9)-Qti{f,9) i^dv 



where h{y) = g{v) + \ Vg{v)\. The kernels and are given by 

B.y{u, u ■ a) = B{u, u ■ (T)|nP, Bj{u, u ■ a) = max(;B(ii, u ■ a), \VuB{u, u ■ a)\)\u\^~^'^ , 
moreover, es{-) is a given restitution belonging to the class Co for any s G [0, 1]. 
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Proof. Using the notation of the Proposition A.|T], we set D^. = le — h- Thanks to ( A..10 ), 
we may spHt De as De = -De,i + ^e,2 with 



De,i = l I f(.v)dv I da I i>{v + z) 
s2 Jn 



dz 



Jei\z\) 



and 



f{v)dv I da I '4'{v + z) 

'§2 JQ* 



Je{\z\) 



Fv,a {'^a{z)) dz. 



We begin estimating |-De,i| which is the more involved part. For fixed v,a, we use the 
following representation formula 

= / VF,^„{U + s{V -U))-Uds 
Jo 

with V = Ce{z) and U = (pcr{z) to get 



Fv,a (Ce(^))-F^,<x i^Az)) 

Therefore 



{fle{z)-l) / (pa{z)-VFy^cr{fa{z) + s{fIe{z)-l)^a{z))ds. 



1 



e,l 



§2 Jn 



f{v)dv / da / i;{v + z) (fieiz) - 1) 



dz 



Je{\z\) 



iPa{z) ■ VF^^a- {^a{z) + s{fj,e{z) - l)(p„{z)) ds 

Now, according to Lemma A. ^, for any s G [0, 1], there exists a restitution coefficient es{-) 
in Co such that fj.es{z) = 1 + s{^e{z) — 1). Therefore, 

iPa{z) + S{^le{z) - l)^a{z) = Ces{z) 

and, performing the backward change of variable z ^ u = C,J^{z) with Jacobian dz = 

Ja{u)Je,{^es ° '^a{u)) du We get 



De,i=4[\s[ \f{v)\dv [ da! ^{v + Ile^o^Au)) 
Jo Jr3 j§2 Jq^ 



/Ue(ng, O $^(m)) - 1 



Since es(-) is a restitution coefficient in the class Co, thanks to the universal bounds ( ]A.90 
we see that 

sup (1 + ti • crj— — — ^ , ,, ^ 16 < oo. 

se(o,i),CTes2 

Moreover, it is easy to see that 



\fJ-eiz) - l| 



Mae{\z\)) 



|1 - I3eiae{\z\))\ ^ i^e) \z\aei\z\y ^ 2^^^(e)|zr+^ Vz G 
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Since llle^ o ^aiu)\ ^ l^'o-l^i)! ^ I""! and \ip„{Ue^ o ^„{u)) \ ^ |n|, we get 



|Z)e,i| < 2^+6^^(e) / \fiv)\dv 



da 



[v., u 



^+^du /V^.,a(n)| ds 
Jo 



where Vs = v + lig^ o ^^(ti). One can check that 

|ii|7+2 |VF„^^(u)| ^ B^{u, u ■ a) {g{v + u) + \Vg{v + u)\) = F^(n, u ■ a)h{v + u). 
Therefore, performing again the change of variable u — > —u, we obtain 



|Z)e,il < 2^+6^^ (e) 



ds 



\f{v)\B^{u, u • a) h{v — u)\il}{v'g)\dv du da 



where v'g = v + o (^^[—u) = v — f3e 



u\ 



is the post-collisional velocity 



associated to the restitution coefficient e^. This proves that 



|De,i| ^ 2^+4^^(e) C ds [ 
Jo Jrs 



Qj , {f,h){v)mv)\dv 



(A. 12) 



where Q- is the collision operator associated to the kernel B-y and the restitution co- 

efficient eg. For the estimate of |De,2| it is enough to prove that there exists Ce > such 
that 

1 



1 



(A.13) 



Indeed, if ( |/V.13p holds then 

\De,2\^Ce [ \f{v)\dv 



da 



{ipiv + z)\ \F^^„ iv^aiz)) I dz. 



Performing the backward change of variables u = ipa{z) as before 



l/Wldt; / da 



{V + «>a('u))| (n) I \u[<J^{u) du 



where we used that |<I>cr(^t)| ^ \u\. Changing again the variable u into —u we get 



e,2 



B^{u,u- a)\f{v)\ \g{Vi,)\ \il^{v[)\ d-udu^ do- 



where v'l is the post-collisional velocity associated to elastic interactions, that is, v[ = v 

u—\u\a p-rii ■ 

— This gives 



\De,2\<Ce Qt,if,g)iv)mv)\dv 

Jm.^ ^' 

which, combined with ( |A.12p yields the result. The idea to prove ( |A.13p is to evaluate 
Jeig) for ~ 0. Since e(r) ~ 1 — ar'^ for r ~ one checks that 



-(l+<(r))/32(r):^l + 



a(7 - 1) 



for 
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Since aeio) ~ for ~ 0, we get Jeio) ~ 1 + '^^^\^q~' as ~ 0. Therefore, 



1 - Je{Q) ^ 

sup < oo 



and ( |A.13P follows for some constant Ce depending only on e(-). □ 

Remark A. 3. Notice that, defining as in Section ^ the rescaled restitution coefficient e\{r) = 
e(Ar) (A > Oj, one sees from the above reasoning that 

Je.(g)^l + ""^^"^^ A^g^ as A^O. 
In particular, for A small enough the constant Ce^ appearing in the above Proposition satisfies 

This property vuill be important in Section I 



Example A. 1. Assume e(-) is the restitution coefficient corresponding to visco-elastic hard- 
spheres 



e(r) = l + r^O. 



k=l 

Then, setting He{r) = ^ (1 + ^e(^)) l^ei'"')' difficult to prove that there is some explicit 

constant C > such that 

\He{r)-l\^C{l-e{r)), Vr ^ 0. 

In particular, \He{r) — 1| ^ C£j{e)r'^ for any r ^ Ofrom Mfhich we deduce that ( \A..1'30 follows 
with a constant Ce proportional to £^{e). Since i^{e\) = X^£j{e), there exists some constant 
c > such that Ce^ ^ cX' for any A G (0, 1] (not just for A small enough as in the previous 
remark). 



A.3. About the energy identity. Recall that, for any solution Gx to ( |1.15p , one has the 
identity 

^Q=T^ f Gx{v)Gx{v,)^e{Av-vA^)dvdv, 

X ^ yM3xR3 

where q = f^,^ G\{v) dv and *e(-) is defined by (]1.9D. Notice that for any fixed r > 0, 

^ *e(AV2) ~ _^r3+T as A ~ 0. 



A3+7 ' 4 + 7 

Define for simplicity 

Ur') = ^*e(AV2) Co(r^) = 4^^'^^ 

and the two functional 

Ixif,9) = / f{v)aiv*)Cx {\v - w^P) dvdv^, 
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and 

We will write Ix{f) = Tx{f, /) and Xo(/) = Xo(/, /). Then, one has the following 

Lemma A. 5. There exist a positive constant vl-y > such that, for any 5 > and e > there 
exists Xq g (0, 1) such that 



sup |Xa(/i,5i) -Xo(/2,52)K -4^ (ll/i - /2||i,i ll^ilLi + lli^i -52||li II/2IIL1 

Ae(0,Ao) ^ 

+ ^(ll/2|Ui|b2|Ui + ||/2|Li^^^Jb2|Li,^J. 

In particular, if g e Ll^^_^_g and / G Ll_^_^, then 

limsup|lA(/)-Xo(3)|^A^||/-3|Li (ll/ll^i +lbllLi )• 

X_^Q 3+7 V 3+7 3+7/ 



Proof Note that \Tx{fi,9i) - Mf2, 92)] ^V^ + Vl + Vl where 

'^1= f 1/1 W - f2{v)\ \gi{v,)\Cx {\v - v,f) dvdv, 
^1= [ I/2WI \91iv.) - 52^)1 Ca {\v - ^*|') dvdv, 

7r3xIR3 

2^A= / |/2(^;)|b2K)| \Cx{\v-V.\^)-Co{\v-V^\'')\dv.dv. 

Let us investigate separately these three terms. Since there is some positive constant 
> such that *e(^^) ^ K^r^+'f, it is clear that Ca {\v - v^l"^) ^ Kj\v - v^^+t ^ 

3+'Y 

2— K^(u)3+->(u^)3+7 for any {v,v^). Therefore 
Vi ^ 2^K^ [ \f,{v) - f2{v)\ \gi{v,)\{vf+^v,f+^dvdv, 

7r3xR3 



2^i^,||/i-/2|Li^Jbi||^i^^. 



In the same way 



Vl^2 — K4g,-g2h^J\f2hl^^. 
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Regarding the term Vf, set ujx{R) = supg^^^i^ ICaC?^^) — CoC?^^)! for any i? > 0. It is clear 
that for any fixed R> one has liniA^o = 0. Let i? > be fixed and spHt as 

-^1= [ \f2{v)\\g2{v.)\\Cx(.\v-v,\^)-(:o{\v-v^\^)\dvdv^ 



+ / \f2{v)\\g2{v,)\ \Cx{\v - v,f) - Co {\v - vX)\ dvdv, 

J\v-v^,\>R 

^Ux{R)\\f2\\L492\\L-^+{K^ + C^) I 1/2^1152^)1 \v - V.f^^ dv dv, 

J\v~v-i,\>R 

where we used the fact that \C\i\v — v^P) — Co (1^^ — 'WitP) | ^ (K-y + Cy)\v — u^p"'''^ for any 
{v, v^). Consequently, for any S > 0, 

Vl ^ UJx{R)\\h\\LA\92\\L^ + I 1/2 Wl \92{V.)\ \V - vA^+'<+' dv dv,, 

R JlR3xM3 

that is, 

Vl^Ux{R)\\f2\\LA\92\W+'r^^^^^\\92\\Ll^^^^^^ 

Taking first i? > large enough and then A small enough we get the conclusion. □ 

Lemma A. 6. Assume that there exist two positive constants a, b > and two exponents 
7 > 7 > such that 

|e(r) — 1 + ar'^'l ^ b r'^ for any r ^ 0. 
Then, there exist two explicit positive constant A^,B^ > such that 

\Ixih,9l)-Uf2,92)\ ^ A (ll/l - f2\\LlJ\9l\\Ll^_^ + \\9l-92\\LlJ\f2\\Ll^^ 

+ B^X-\\f2hi ^_||92||li ^_ VA G (0,1) (A.14) 

3+7+7 3+7+7 

where a = min(7,7 — 7). 

Proof. For any A G (0, 1] and r > 

Then, under our assumption on e(-), there are three constants A, i?, C > such that 

\C\{r'^) - Co(r^)| ^ A\~-^r^+^ + BX^ r^+'^^ + C7A~r3+T+" VA > 0, r > 0. 
In other words, there is > such that 



Cx{\v-vA'')-Qo{\v-vA'') 



where a = min(7, 7 — 7). Consequently, 

Vl^C,\^f2\\Li^ \\g2\\Ll^ . (A.15) 

3+7+7 3+7+7 

With this estimate the proof follows as the proof of Lemma A.^ □ 
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Remark 4.11. For visco-elastic hard-spheres, the assumption (4^9) is met with 7 = and 
7 = |. 7n particular, a = ^• 

For a given a ^ define the exponential weight ma{v) = exp(a|ti|) and introduce 

X = L\ma) and y = L\{ma). 

Define 

X = Xr\\f: [ f{v)dv= [ f{v)vdv = o\ and ^ = ^n|/ : / f {v)\v\'^ dv = o\ , 
[ Jr3 Jm_3 J I Jr3 J 

and the operator 

A : h£ X ^ Ah = {Aih;A2h) £Rx X 

where 

Aih = 2Io{M,h) and A2h = Qi{h,M) + Qi{M,h) = ^ih. 
The operator ^ is a suitable lifting operator of . 

Lemma A. 7. The linear functional 

A -.y — >RxX 

is invertible and the norm ||^~^|| = ll-^'^llKx^-i-j) estimated explicitly. 

Proof. The fact that the mapping A2 = .^i ■ X — > X is invertible with explicit inverse is 



a direct consequence of Proposition [4.3| (see []2T|] for details). Set 



p^:= [ {\v,\^ -3e)M{v)M{vM\v-v*\^)<ivdv* 

= C^ [ {Iv^l"^ - Se) M{v)M{v^)\v - v^f'+^ dvdv^. 



Direct inspection shows that / for any 7 > 0. Arguing as in Lemma 4.3], we 
deduce that, for any y G M, 51 G the unique solution to the equation Ah = {y, g) is given 
hy h = hi (fi + with 

= iff hi = -^(y- Aih^) . 

This proves the Lemma. □ 



Appendix B: Existence of a steady solution for diffusively driven granular gases 



The main objective of this section is to prove Theorem |1.2| , that is, to prove the exis- 
tence of an steady solution F to Ql.lp . The proof, see , is based on a dynamic version 
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of Tykhonov fixed point theorem and it is achieved by controlling the L^-norm, the mo- 
ments and the regularity of the solution to the time-dependent problem associated to ( Jl . Ip . 
Consider the diffusively driven Boltzmann equation 

dtfit, v) = Qeif, m, v) + f,Af{v, t) t>0,veR^ 
f{0,v) = fo{v) veW", (B.l) 

with /i > and where the initial datum /o is a nonnegative velocity distribution satisfying 

/ fo{v)dv = l, [ fo{v)vdv = and / fo{v)\vf'dv<oo. (B.2) 

Notice that if Ef{t) denotes the kinetic energy of f{t,v) at time t ^ 0, that is, Ef{t) = 
/k3 fit,v)\v\'^dv then it satisfies 



dt 



Ef{t) = -Uf{t))+6^i 



where le is the energy dissipation functional defined by ( I.IC ) (justifying, a posteriori, the 
terminology we used in the core of the paper). Problem ( B.ip is well posed due to the 
following theorem. 



Theorem B. 1. Assume the restitution coejficient e{-) satisfies Assumption \1 . 1\ and the initial 
datum fo e L^(M'^) n LlogL(M'^) satisfies ( |B.2p . Then, there exists a unique nonnegative 
weak solution 

f G L~([0,oo),L1(IR3)), /log/ G L~([0,oo),L1(M3)) 



to equation ( |B.ip , with the initial condition f{-,0) = fo- Furthermore, if in addition fo G 
thenf£ C^{[to,oo),S{R^))for every to > 0. 



LlnL^ 



The proof of Theorem ^]2| can be deduced from Theorem B.|l| following the proof of 
[Q, Theorem 5.2], thus, we shall only recall the main steps in the proof of Theorem B.|l[ 
The proof will follow the path presented in [ ]T4l , Theorem 5.1] with the differences clearly 
explained. 



B.l Povzner-type inequalities. We derived in m\ Povzner's estimates in the spirit of [11] 



and [19]. We shall extend this result, using some ideas of [14]. Recall that for any non- 
negative function / and text function ^lJ{v) = ^dvl"^) with ^ nondecreasing and convex 



fiv)fiv.)<l>i\u\)AB,emiv^V*) dv.dv, 



where 



= A+jnv,v.)-{^i\v\') + ^{\vX)). 
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The collision cross-section B{u,a) is given by ( |A.lp with the normalization assumption 
( |A.3| ) (notation is slightly changed with respect to [^). Define the velocity of the center of 

mass U = so that 



\u\ 



U 



u\ 



U + Y^ and < . 2 



with w = (1 — /3){i + f3a. We proved in Eq. (2.15)] that the post-coUisional integral can 
be estimated from above as follows 



{C/cr>0} 



^ E 



3 + U-a 



l-U -a 



b{u ■ a) da. 



Under assumption ( |A.4| ) one can prove, see [11, Lemma 1], that this integral (involving tj 
and u) takes its maximum value whenever U = u, that is, 

-1 



^ E 



,3 + s 



.1 



b{s) ds 



(B.3) 



where b{s) = b{s) + b{—s) and E = + jw^p. At this point, we shall adopt the viewpoint 
of [14] and assume that ^ satisfies the following conditions: 

^{x) ^0, X > 0; ^(0) = 0; (B.4a) 
^! is convex, ^" G Lf^^HO, oo)); (B.4b) 
^''(ax) ^ r]i{a) ^''(x) and ^"(ax) ^ 772(a) ^"(x), x > a > 1, (B.4c) 

where 7?i(-) and 7/2(-) are locally bounded functions. Then, one has the following general- 
ization of []T4l, Lemma 3.3] to non constant restitution coefficient. 



Proposition B. 1. Assume that ^'(x) satisfies (B.4). Then, for any {v,v^) G M'^ x M , 

where A = rji{2) while k > is an expUcit constant depending only on r/2 and on b{-). For 
instance in the hard-sphere case b{-) = then k 7?2(2) = 

Proof. Recall that, see [Q Lemma 3.1], if ^ satisfies ( pTl ) then 

^'(x + y)- ^'(x) - -^{y) ^ ^ ( X ^''(y) + y ^''(x) ) (B.5) 

and 

*(x + y) - *(x) - ^'(y) ^ ao xy ^"(x + y), (B.6) 

where A = 771(2) and ao = (2r/2(2))^^. Let (f,w^) be fixed and write AB,e[^]{v.,Vi,) = 
'P[^]{v,Vi,) -7\A[^'](D,ti^) where 



Vm{v, V,) = j^^ l^{\v\^ + \v.f) - *(|H') - *(|^.|')) 



biu ■ a) da 



and 



^{\v\^ + \v^\') - ^{\v'\'^) - ^(Kl^) ) b{u ■ a) da. 
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Using (B^) and the normalization assumption ( |A.3[ ) one gets directly that 

Vm{v,v,) ^ A {\v\H'{\vX) + \vX^'{\v\^)) . 
Let us extimate A/'[^'](f , v^) from below. First, one notices that 

and deduces from ( |B.3p that 



Second, since b{s) ds = ^ according to ( lA.3p one can write 



b{s)ds, S=|up + |v^p. 



b{s)ds. 



Noticing that E = E^^ + E^-j^ for any s £ (0, 1), it is possible to apply directly (pjH) to 
obtain 

Af[^]{v, v^) ^ ^eH"{E) / (3 + - s)b{s) ds. 
Setting k = ^ (3 + — s)b{s) ds, the desired conclusion follows. 



□ 



Remark B. 1. Note that the above estimate does not depend on the restitution coefficient e(-). 
Indeed, the two constants A and b are depending only on ^ and the angular cross-section b{-) 
but not on e(-). 



With Proposition B.|l| the following properties are derived exactly as shown in []T4|]. 
Lemma B. 1. Let p > 1 and ^(x) = x^. Then, for b{-) = I/Att one has 



\v-V,\lCem{v,V.) ^ -kp{\v\^P+^ + 



|2p+l 



where the constants kp and Ap are independent on the restitution coefficient e(-). As a conse- 
quence, for any nonnegative distribution f = f{v) ^ 0, 

Qe(/, f){v) \v\''' dv ^ -kp f{v) dv^ f{v)\v\^P+' dv^ 

fiv)\v\dv) ( [ fiv)\v\^Pdv 



B.2 Proof of Theorem B. |l|. The proof is a modification of [Q, Theorem 5.2] and we only 
give a sketch of it explaining where the original argument has to be modified to handle the 
non-constant restitution coefficient. Using our Povzner's estimates, the propagation and 



appearance of moments given in []Tj, Lemma 3.5] follow. Additionally, using the control 
of \\QB,e{f, DWlp derived in Corollary A. we can easily adapt the proof of [|T^ Lemma 
4.7] to our case, yielding a local in time propagation of H^(]R^) norms. Therefore, the a 



priori estimates for the solution to (B.l) derived for the constant restitution case extends. 
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Let us first deal with a smooth initial datum /o with compact support. For any truncation 
parameters M > 1 > m > 0, define then 

^m,M(|u|) = m + min(|ti|,M) 

and set Bm,M{u,a) = ■^^m,M{\u\), u G M^, a G S^. Define the collision opera- 



tor Qm,M = QBmM,e (using the notations of equation ( lA.2p ). For any T > 0, let 



g = g(t,v) G L°°([0,T] ; Ll{M.^) D L'^{M.^)) be a nonnegative function with 

/ g{t,v)dv = l and / g{t,v)dv = £ [0,T]. 

Jr3 Jr3 

Consider the auxiliary problem 

f dtfit,v) - fiAJit,v) + Mf{t,v) = Qm,Mi9,9)it,v) + Mg{t,v) t e [0,T] , v e M.^ 
\ f{0,v) = Uv). 

(B.7) 



Setting h = Qm,Mig,9)it,v) + Mg{t,v), one checks, see [14, Theorem 5.2], that h £ 



L°°{[0,T] ; L^(M3) n L2(m3)) and h ^ -g{g * ^>m,A/) + Mg ^ 0. The unique solution 



/ G L°°([0, T] ; L^(M3)nL2(M3)) ( |bj| ) can be given explicitly and by a classical parabolic 
regularity result 

II/IIh2{[o,t]xr3) ^ Cm(||^|| R3) + II/o||hi(R3))- (B.8) 

Denoting by T the operator that maps g into /, the core of the proof consists in showing 
that for a certain choice of constants Ai and A2, the operator T maps B into itself. Here 
we refer to the set, 

B= \f eL\[0,T]xR^) ■ f^O, [ f{t,v)dv = l, [ fit,v)vdv = 0, 

I JR3 g-j 

Ef{t):=[ f{t,v)\v\'^dv^Ai, [ f{t,v)dv^A2 for a.e. t G [0, T] I . 

The first three properties are clearly satisfied. To determine Ai, one multiplies equation 
( |B.7[ ) by |f p and integrate by parts. This yields 

d f 

—Ef{t) + MEf{t)^6fi + MEg{t)+ / Q™,M(5,5)(i,^^)kl^dt; 

(^6fi + M^Eg{t)- / g(t,z;)g(t,z;,)^^,M(|-»|) , ^ dt;dz;,, 

yR3xM3 \U\ 



where we used ( ]1.8[ ) in the last identity. Since ^dr) may be arbitrarily small for small 
r > 0, the argument changes slightly with respect to [^^. Using that limsup^^oo e(^) = 
eo < 1, there exists some iio ^ 1 and some constant C > such that 

^e{\u\'^)>C\uf \/\u\>Ro. 
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Therefore, 

g{t,v)g{t,v^)^rnMiH)—^TT~^ dwdu* ^ C [ g{t,v)g{t,v^)^rn,MiH)W\'^ dvdv^ 

'xR3 \U\ J\u\^Ro 

^ Cm g{t, v)g{t, v^:)\u\'^ dvdv^. 

J\u\^Ro 

Since g has unit mass, 

I git, v)g{t, v^,)\u\'^ dvdv^: = / g{t,v)g{t,v^:)\u\'^ dv dv^, 

- / g{t,v)g{t,v^)\u\'^ dvdv^ ^ 2Eg{t) - Rl 

J\u\<R.o 

Going back to ( |B.10p finally leads to the estimate 

-^Ef{t) + MEf{t) ^ 6^ + MEg{t) - 2CmEg{t) + CmRl (B.ll) 

Setting A[ = + CmRQ assuming Eg{t) ^ A'^ 3delds the differential inequality 

-^^Ef{t) + MEf{t) ^ MA[ 
which, in turn, implies Ef{t) ^ max {A[,Ef{0)) for any t e [0, T]. Thus, one may choose 

Ai = max (^i, foiv)\v\^ dv^ 

in the definition ( |B.9P of B. For the determination of the parameter ^2 > just follow 
the path of [ |l4i Theorem 5.2]. This leads to the existence of a solution / = fm,M £ 
L°°([0, T],Ll{R^) n L2(m3)) to the modified Boltzmann equation 

idtf{t,v) =Qm,MifJ)ii^^) + t^^vfit,v) t£[0,T],vGR^ 
\f{0,v) =Mv). 

It remains to pass to the limit as M — )- oo and m — ;> 0. To this end, we will show that the 
bounds found in the a priori estimates hold for the fixed point solutions and are uniform 
in M and m. From ( [B.lip with / = fm,M 

-^Ef{t) ^ 6/z - 2CmEf{t) + CmRl, 

which yields 

which provides a bound independent of M. Using Proposition Bjl], it is possible to adapt 
the proof of [Q, Theorem 5.2] to get that, for any p ^ I and T > 0, the bounds of 
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/ = fm,M in L°°([0,r],L|p(M3)) are independent of M. Since / e H2([0,r] x M^), using 
the extension of Lemma 4.7] and then [jT4|, Lemmas 4.8 & 4.9], 

feL^{[0,T],Ul{R')), 

for every n ^ 1, and every p ^ 0, with bounds independent on M. This allows to pass 
to the limit as M — ^ oo in the weak form and to show that the limit solutions satisfy the 
equation with the kernel 

(m + \u\) b{u, a). 

Following the argument of [jl^] it is possible to prove that the bounds in 
L°°{[0,T],Llp{M.^)) are actually independent on m and T. This allows to pass to the limit 
as m — > and the limit solution obtained is a solution to ( jB.lp . A standard approximation 
argument generalize the initial conditions from smooth ones. 
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